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PREFACE 

r  I  "HE  Theory  of  Map  Projections  has  received  the  attention  of 
-•-  so  many  eminent  mathematicians  and  geodesists  that  it  may 
seem  somewhat  presumptuous  to  offer  anything  further  on  the 
subject.  These  Papers  must  therefore  prove  their  own  justification. 
They  are  the  result  of  investigations  extending  over  (with  intervals) 
nearly  four  years,  by  one  who  was  for  twenty  years  employed  on 
the  Trigonometrical  and  General  Survey  of  one  of  the  British  Pro- 
tectorates in  Tropical  Asia,  but  who  had  until  recently  devoted  little 
attention  to  the  theory  of  map  projections,  being  content  to  follow 
the  established  methods  of  the  Ordnance  Survey  and  the  Survey  of 
India. 

The  investigations  came  about  in  this  way.  At  the  beginning 
of  1916  the  writer  had  offered  his  services  to  the  Secretary  of  the 
Royal  Geographical  Society  to  assist  in  the  map  work  being  done 
by  that  Society  for  the  War  Office.  One  of  the  first  tasks  given  him 
by  Mr  Hinks  was  to  compute  some  Tables  of  rectangular  coordinates 
of  the  special  modified  Polyconic  Projection  adopted  for  the  Inter- 
national  Map  on  the  scale  of  I/million,  for  a  few  sheets  to  be 
compiled  on  half  that  scale.  The  carrying  out  of  this  calculation 
led  to  an  examination  of  the  basis  of  the  formulae  on  which  the 
modification  referred  to  depends.  This  led  to  a  comparison  with 
other  Projections  and  their  errors,  until,  becoming  completely 
fascinated  with  the  subject,  the  writer  thinks  he  can  claim  to  have 
passed  under  review  the  principal  books  and  Papers  which  have 
been  published  on  it,  and  also  to  have  worked  out  some  results 
which  are  new  so  far  as  he  knows,  and  which  he  hopes  may  prove 
useful  to  practical  map-makers  and  surveyors.  These  results  are 
concerned  mostly  with  the  Minimum  Error  Projections,  and  although 
the  writer  cannot  claim  to  have  discovered  anything  of  striking 
novelty  in  either  theory  or  practice,  he  has  been  able  to  extend 
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the  method  invented  by  Airy  and  Clarke  more  generally,  and  to 
simplify  their  results  ;  so  that,  out  of  the  almost  bewildering  number 
of  Projections  that  have  been  discovered  and  advocated  from  time 
to  time,  those  which  are  practically  the  best  or  most  useful  are 
reduced  to  comparatively  few,  the  application  of  which  has  been 
simplified.  One  of  the  most  interesting  results  is  to  establish  in 
their  proper  place  some  old  and,  by  later  writers,  somewhat  dis- 
paraged Projections,  which  seem  to  have  been  hit  upon  by  their 
inventors  by  a  sort  of  geometrical  intuition  rather  than  by  rigorous 
analysis. 

Owing  to  the  present  high  cost  of  printing  the  Papers  are  pre- 
sented in  as  concise  a  form  as  possible.  They  do  not  profess  to  be 
a  complete  treatment  of  the  subject  ab  initio,  but  only  a  statement 
of  what  are  supposed  to  be  new  results.  Many  steps  in  the  mathe- 
matical processes  have  been  left  out,  and  in  the  last  parts  it  will  be 
difficult  to  follow  the  argument  without  reference  to  the  works  of 
the  other  writers  mentioned. 

The  Tables  etc.  have  been  calculated  with  5 -figure  logarithms 
and  Fuller's  slide-rule,  and  are  intended  only  to  be  illustrative.  In 
the  formulae  the  symbol  log  means  the  natural  logarithm  and  the 
angles  are  supposed  to  be  reckoned  in  circular  measure.  Though 
the  work  has  been  checked  where  possible  by  alternative  methods, 
there  has  been  no  absolutely  independent  check,  and  it  cannot  be 
that  some  errors  have  not  crept  in  both  in  formulae  and  arithmetical 
work.  The  writer  will  therefore  welcome  corrections  from  any  reader 
who  may  have  the  will  and  opportunity  to  check  his  results. 

In  the  course  of  these  investigations  the  writer  received  valuable 
hints  as  to  suitable  and  fruitful  lines  of  inquiry  from  Mr  A.  R.  1  links, 
F.R.S.,  Secretary,  Royal  Geographical  Society,  to  whose  own  little 
book  on  Map  Projections  he  is  indebted  for  a  clear  preliminary 
survey  of  the  whole  subject.  He  has  also  received  valuable  assistance 
in  mathematical  difficulties  from  Professor  A.  F.  II.  Love,  F.R.S., 
Oxford,  Dr  T.  J.  I'A.  Bromwich,  F.R.S.,  Cambridge,  and  Professor 
J.  \Y.  Nicholson,  F.R.S.,  King's  College,  London.  lie  is  indebted 
to  Captain  G.  T.  McCaw,  General  Staff,  for  reading  the  original 
draft  and  for  making  some  useful  suggestions. 


PREFACE  vii 

The  following  is  a  list  of  the  principal  works  consulted.  When 
quoted  in  the  text  they  will  be  referred  to  by  author's  name  and 
the  page,  e.g.  Germain,  p.  150,  for  the  sake  of  brevity  : 

GERMAIN,  A.    Traite  des  Projections  des  Cartes  Geographiques,  Paris,  1865. 
FIORINI,  M.    Le  Projezioni  delle  Carte  Geografiche,  Bologna,  1881. 
TISSOT,  A.    Memoire  sur  la  Representations  des  Surfaces  et  les  Projections 

des  Cartes  Geographiques,  Paris,  1881. 

CRAIG,  THOMAS.    A  Treatise  on  Map  Projections,  Washington,  1882. 
HERZ,  NORBERT.    Lehrbuch  der  Landkartenprojektionen,  Vienna,  1885. 
HELMERT,  E.  R.    Hohere  Geodasie,  Vol.  i,  Leipzig,  1880. 
HAMMER,   E.     Uber  die  geographischen  wichtigsten  Kartenprojektionen, 

Stuttgart,  1889. 

PIZZETTI,  P.    Trattato  di  Geodesia  Teoretica,  Bologna,  1905. 
CRAIG,  J.  I.    The  Theory  of  Map  Projections,  Cairo,  1910. 
HINKS,  A.  R.    Map  Projections,  Cambridge,  1912. 
ADAMS,  OSCAR  S.   Special  Publications  of  United  States  Coast  and  Geodetic 

Survey,  Nos.  47,  49,  52,  57,  Washington,  1917-19. 
Encyclopaedia Britannica.  Articles  on  Mathematical  Geography,  gth  Edition, 

and  on  Map  Projections,  nth  Edition. 
JORDAN,  W.     Handbuch  der   Vermessungskunde,  Vol.   in,    6th    edition, 

Stuttgart,  1916. 

A.  E.  YOUNG. 

ROYAL  GEOGRAPHICAL  SOCIETY, 

KENSINGTON  CORK,  LONDON,  S.W.  7. 
May  1920. 


CONTENTS 

PAGE 

Chap.  I.    The  Minimum  Error  Zenithal  Projections         .        .  I 

Airy's  Zenithal  Projection    ........  2 

Breusing's  Projection  .........  7 

Orthornorphic  and  Equal  Area  ....... 

Equidistant           ..........  9 

With  total  area  true     .........  10 

Perspective           ..........  12 

With  scale  true  at  centre    .         .         .         .         .         .         .         .  17 

Evaluation  of  Average  Error       .......  19 

Chap.  II.    The  Minimum  Error  Conical  Projections        .        .  22 

Orthornorphic      ..........  23 

Approximately  Equal  Area          .         .         .         .         .         .         .  25 

Minimum  Error  Map  of  a  Zone           ......  26 

Radial  scale  true  ......... 

Plate-Carree         .......... 

Cassini          ........... 

Merrator 

Conformal  of  Gauss     .........  35 

Euler's  with  two  standard  parallels     ......  35 

Murdoch's  Projections         ........  38 

Equidistant  with  total  area  true           .         .         .         .         .         .  41 

Albers'  Equal  Area      .........  4; 

Collision's  World  Map       ........  48 

G.  W.  Hill's  Projection        ........  40 

Pjest  projection  for  given  country        ......  49 

Mathematical  paradox         .         .         .         .         .         .         .         .  - - 

Chap.  III.    The  Conical  Orthomorphic  Projection  with  two 

standard  parallels  (Lambert's  second)  for  the  Spheroid.  57 

Chap.  IV.    The  Polyconic  Projections   ....  (^ 


Chap.  V.    The  finite  errors  of  Projections     . 

Chap.  VI.    '1  he  Conver^cncy  of  Meridians  . 

Expansions  used  in  the  work 


CHAPTER  I 

THE  MINIMUM  ERROR  ZENITHAL  PROJECTIONS 

The  Zenithal  Projections  are  a  particular  case  of  the  Conical  Class 
obtained  by  opening  out  the  cone  until  it  becomes  a  plane  or,  which  is  the 
same  thing,  by  making  the  constant  of  the  cone  equal  to  unity.  They  will 
be  considered  before  the  general  conical  class,  however,  because  they  are 
the  first  to  which  the  principle  of  Minimum  Error  was  applied,  and  it  will 
be  simpler  to  approach  the  subject  from  what  has  already  been  done. 

Suppose  that  r  is  the  distance  from  the  centre  or  pole  to  any  point  on 
the  map,  then  if  0  is  the  corresponding  angle  subtended  by  this  distance 
at  the  centre  of  the  Earth,  supposed  spherical,  we  have  as  the  law  of  the 
projection  r=JRf(6\  where  R  is  the  radius  of  the  Earth.  In  the  following 
investigations  we  shall  for  simplicity  assume  R-  unity,  but  it  must  be 
remembered  that  to  convert  the  values  of  r  in  any  case  into  actual  dis- 
tances, they  must  be  multiplied  by  the  proper  value  of  R  on  the  appropriate 
scale. 

The  radial  scale  of  the  map  is  given  at  once  by  drfdO  and  the  perpen- 
dicular scale  by  r/s'm  6. 

By  choosing  different  forms  of  the  function  r—f(&)  we  obtain  Zenithal 
Projections  with  different  properties  of  which  the  following  are  the  principal : 

Name  r=f(ff)  Radial  Scale  Perpendicular  Scale 

=  dr\dQ  =r/sin0 

1.  Equidistant  0  i  cosec  6 

2.  Orthomorphic 

or  Stereographic          2  tan  J0 

3.  Equal  Area  2  sin  |-0 

4.  Gnomonic  tan  0 

ksmO 
^.    Perspective  , 

h  +  cos  0 

It  will  be  observed  that  the  fifth  only  is  multiplied  by  a  constant  /-, 
commonly  called  a  scale  constant,  but  it  will  be  shewn  later  that  the  proper 
evaluation  of  this  constant  plays  a  most  important  part  in  the  determination 
of  the  criterion  of  minimum  error,  and  that  a  similar  constant  could  be 
applied  to  the  other  four  with  like  result,  which,  so  far  as  the  writer  is  aware, 
has  not  hitherto  been  noticed.  The  second  projection  in  the  table  has  the 
property  that  the  scale  in  any  direction  round  a  point  is  the  same,  and  the 
third  the  property  that  the  radial  and  perpendicular  scales  at  any  point  are 
reciprocal,  so  that  equality  of  small  areas  round  any  point  is  preserved. 
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Airy's  Zenithal  Projection  by  Balance  of  Errors. 

Sir  George  Airy,  late  Astronomer  Royal,  conceived  the  idea  of  a  pro- 
jection which  should  be  a  sort  of  compromise  between  these  two.  His 
investigation  is  published  in  the  Philosophical  Magazine  for  December 
1861.  It  consists  in  finding  by  mathematical  analysis  the  form  of  the  func- 
tion r  --•/(&)  which  makes  the  squares  of  the  two  scale  errors,  integrated 
over  the  surface  of  the  sphere  represented  by  the  map,  a  minimum.  Airy's 
method  consists  in  fact  of  making  the  average  scale  error,  or  more  strictly 
speaking  the  mean  square  scale  error,  of  the  map  a  minimum  ;  this  average 
being  found  by  squaring  the  scale  errors  in  the  two  directions  mutually 
perpendicular,  integrating  them  over  that  portion  of  the  sphere  represented 
by  the  map,  then  (implicitly)  dividing  by  the  area  of  the  same  surface  and 
extracting  the  square  root.  The  result  will  be  the  mean  square  scale  error 
of  the  map,  from  which  the  average  scale  error  can  be  found  by  multiplying 
by  07979.... 

Unless  orthomorphisni  or  equality  of  small  areas  is  wanted  for  some 
particular  purpose,  it  is  difficult  to  conceive  a  more  rational  criterion  than 
Airy's  for  finding  the  best  projection.  Another  criterion  is  to  make  the 
maximum  positive  scale  error  equal  to  the  maximum  negative  scale  error, 
in  other  words  to  make  the  maximum  scale  error  as  small  as  possible. 
Since  the  scale  error  increases  as  the  square  of  the  distance  from  the  pole 
of  projection,  which  is  the  centre  of  the  map  in  the  zenithal  class,  this  has 
the  effect  of  diminishing  the  scale  error  at  the  margin  of  the  map  at  the 
expense  of  increasing  it  near  the  centre  and  increasing  therefore  the  average. 

Nevertheless,  Airy's  Projection  by  Balance  of  Errors  has  not  been  much 
used,  and  it  seems  to  have  attracted  more  attention  on  the  Continent  than 
in  this  country.  The  reason,  probably,  is  the  apparent  complexity  of  its 
formula,  but  also,  the  writer  thinks,  because  its  nature  has  not  been  fully 
understood  or  explained  in  the  text  books — the  more  so  since  Airy  himself 
made  a  slip  in  his  solution  which  caused  the  projection  to  appear  under  a 
disadvantage  in  the  Tables  published  with  it.  Although  this  slip  was  dis- 
covered by  Clarke  and  rectified  in  a  paper  published  in  the  Philosophical 
Magazine  for  April  1862,  Airy's  erroneous  tables  have  unfortunately  been 
incorporated  in  later  publications,  e.g.  A  Treatise  on  Projections  by  Thomas 
Craig,  United  States  Coast  and  Geodetic  Survey  (Washington,  1882),  and 
others. 

It  also  seems  to  be  implied  by  some  writers  that  Airy's  is  the  Minimum 
Error  Projection,  whereas  it  is  only  the  Minimum  Error  Projection  of  the 
Zenithal  Class,  which  is  a  particular  case  of  the  more  general  Conical  ('lass. 
The  writer  lias  succeeded  in  applying  Airy's  method  to  this  more  general 
class,  and  the  results  derived  will  be  given  later. 

( iermain  gives  a  full  account  of  Airy's  Projection  with  Clarke's  correction 
and  in  his  general  review  p.  22"]  says,  "  La  projection  de  sir  Airy,  Halunce 
of  Errors,  n'a  pas  encore  ete  employee;  sa  formule,  tres-complique'e  neces- 
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site  des  calculs  tres-penibles,  et  les  problemes  de  la  sphere  ne  peuvent 
s'y  resoudre  facilement."  Its  qualities  and  the  modification  of  it  called 
Breusing's  Projection  are  also  treated  by  Hammer;  but  the  writer  has 
found  by  far  the  most  complete  treatment,  with  some  extensions  of  Airy's 
and  Clarke's  method,  in'Le  Projezioni  delle  Carte  Geografiche  by  Pro- 
fessor Matteo  Fiorini  (Bologna,  1881),  for  reference  to  which  he  is  indebted 
to  a  pamphlet  by  J.  I.  Craig  of  the  Egyptian  Survey  Department  on  The 
Theory  of  Map  Projections  (Cairo,  1910). 

Returning   to    Airy's    solution,    the   two   scale   errors   at   any   zone   at 
angular  distance  0  from  the  centre  of  map  are  drjdO-i  and  r/sinti—i. 
The  corresponding  element  of  area  on  the  sphere  is 
RdB  .  2irR  sin  6  =  2-irR-  sin  8d6, 

and  the  problem  therefore  is,  omitting  the  constant  multiplier  277- /f2,  to  find 
the  function  r--f(&)  which  will  make  the  integral 

p   (fdr       \-      /    r          \~\ 

M  —  \     \  I  -,„  —  i  }  +  [  — — -,—  i  )  }-  sin  0  dQ  a  minimum, 
Jo  \\dB      J      \sin0      /  J 

where  ft  is  the  angular  radius  at  the  limit  of  the  map.    The  corresponding 

;'£ 
area  on.  the  sphere  is  2irR-  I    sin  Q  d9  =  i  —  cos  ft,  omitting  the  constant 

multiplier  zirR- ;  so  that  the  square  of  the  mean  square  error,  which  we 
shall  call  /u.2,  is  given  by  Mj{i  -~ cos/3).  By  the  Calculus  of  Variations,  Airy 
finds 

r  =  A  tan  |-0  +  2  cot  |  (9  log  sec  W  +  B  cot  \Q, 

where  A  and  B  are  constants  to  be  determined  by  the  boundary  conditions. 
Since  r  must  =  o  when  6  =  o,  Airy  and  Clarke  put  at  once  B  =  Q.  Airy  then 
makes  the  mistake  of  supposing  that  the  constant  A  is  completely  at  his 
disposal  and  he  chooses  A  so  that  the  scale  error  shall  vanish  at  the  centre, 
i.e.  he  takes  A  —  i,  so  that 

r  -  tan  \Q  +  2  cot  \6  log  sec  J#  simply, 

from  which  his  Tables  have  been  calculated.  The  Calculus  of  Variations 
exacts  however  that  at  each  boundary  drjdO-i,  and  we  ought,  strictly 
speaking,  to  suppose  the  map  to  extend  from  a  lower  limit  B-a.  to  a  higher 
limit  0  =  ft,  at  each  of  which  drjdB=  i. 

From  these  two  equations  we  should  obtain 
_     /log  sec  lift  —  log  sec  |a\ 

2  \      tan-  |y8  -  tan-  \a     J  ' 
_     /tan2  |a  log  sec  |-/2  -  tan2  |/?Jog  secj_a\ 
2  \  tan2 1/3- tan2  ^ a 

If  then  our  map  extends  from  the  centre  we  put  a  =  o  and  find 

A-  2  cot2  i/3  log  sec  -\ft,     B  =  o, 

which  are  the  values  of  the  constants  found  by  Clarke.  If  our  map  were  to 
extend  over  a  zone  of  the  Earth's  surface  from  a  spherical  radius  a  to  a 
spherical  radius  ft,  we  should  reduce  the  average  scale  error  by  making  use 

i "> 


4  MINIMUM  ERROR 

of  the  constant  B  as  calculated  above,  and  the  term  B  cot  \Q.  It  is  not 
likely,  however,  that  such  a  map  would  be  of  much  use  in  the  zenithal 
system,  though,  as  we  shall  see  later,  this  term  is  of  great  importance  in  the 
zonal  maps  of  the  Conical  Class. 

Let  us  now  examine  the  separate  terms  of  Airy's  solution.  The  first, 
A  tan  i#,  is  evidently  half  the  characteristic  formula  of  the  Orthomorphic, 
r  =  2  tan  1>6,  or  it  is  the  solution  of  the  differential  equation  drfdO  —  r/sm  0. 
The  second  appears  more  complex,  but  if  we  differentiate  it  we  find 

—  (2  cot  16  log  sec  \0)  -  i  —  cosec'-'  i$  log  sec  \6, 

and  if  we  divide  it  by  sin  6  we  find 

2  cot  \6  log  sec  |$/sin  6  =  cosec2  \B  log  sec  -^6, 

so  that  2  cot  \Q  log  sec  \B  is  half  that  function  which  has  the  property  that 
its  radial  scale  error  is  equal  and  opposite  to  its  perpendicular  scale  error. 
It  is  in  fact  half  the  function  which  solves  the  differential  equation 

(dr  d&  -  i)  +  (r/ sin  6  -  i)  ^  o, 

that  is,  r  —  4  cot  \B  log  sec  \Q  +  B  cot  \  6, 

where  in  the  present  case  B  —  o.  It  is  therefore  to  a  first  approximation 
half  a  function  giving  an  Equal  Area  Projection,  the  characteristic  function 
of  which  is  r=  2  sin  \B.  This  is  seen  more  readily  from  the  expansions  of 
these  two  functions,  viz., 

aani0  =  0-^» +  „*„*" -»•, 
4  cot  W  log  sec  J0  =  <0  -  ,'4.0»  -  Tio  &  -•••• 

We  see  therefore  that  Airy's  Projection  is  in  form  practically  a  mean 
between  an  Orthomorphic  and  an  Equal  Area  Projection,  but  the  ortho- 
morphic  part  must  be  multiplied  by  a  constant  A  which  has  the  effect  of 
reducing  the  positive  scale  error  at  the  outer  boundary,  at  the  expense  of 
introducing  a  negative  scale  error  at  the  centre.  The  function 

r—  4  cot  \B  log  sec  \Q 

does  not  seem  to  have  been  used  by  itself  for  a  map  projection,  nor  does 
there  seem  to  be  any  reason  for  preferring  it  to  the  simpler  and  practically 
identical  Equal  Area  function  r-2sini$.  It  lends  itself  more  readily, 
however,  to  the  mathematical  operations  necessary  to  determine  the  con- 
stants of  the  Minimum  Error  Projection  and  is  a  useful  guide  when  dealing 
with  the  absolute  Equal  Area  function  in  both  the  Zenithal  and  Conical 
Classes. 

By  substituting  the  value  of  r  as  found  above  in  the  expression  for  M 
and  integrating,  we  can  find  the  value  of  M  and  hence  of  //.  in  terms  of  ft. 

We  shall  find 

S  (log  sec  1/J)' 
M     zsin8^-          t;^^         ' 

4  (log  sec  .',£)-  sec  \3 
,¥  (i -cos/3):    ,_  f- 
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Clarke  in  his  Paper  quoted  above  does  not  appear  to  have  taken  this  step, 
but  he  evaluates  the  value  of  M  for  the  particular  case  of  /3=  113°  30'  by 
mechanical  integration,  getting  the  value  M=  0-1569.  The  above  formula 
gives  0-15688,  which  checks  Clarke's  value  exactly. 

Airy's  Solution  expanded. 

It  is  instructive  and  more  useful  for  the  purpose  of  comparison  with 
other  solutions  to  expand  Airy's  solution  as  follows  : 

We  have  r  ~  A  tan     6+2  cot  -  B  log  sec     B. 

222 

Now  tan  -  6  =  -  B  +  —  O3  +  —  05  +  — 7—  B7  +   . 

2  2  24  240  40320 

2  cot  -  0  log  sec  -  0  -  -  6  -  \  B*  -    \    Br> l—  B'  - 

2  2  2  48  960  17920 

•f— 06 —    (f  +    .... 

I  2O  504O 

We  know  that  the  first  term  of  the  constant  A  expanded  in  terms  of  (3  is  i, 
and  we  shall  save  some  trouble  by  writing  A  =  i  + 1>,  where  b  only  is  to  be 
found.  We  can  then  add  the  first  two  expansions  together  and  write 

i  i  59/j-  * 

48         320         161,280  2 

thus 

4  24 

i  i 

4  24 
therefore 


M=  f  *  {(drldB  -  i)-  +  (r/sin  0  -  if}  si 
Jo 


JO    (256       V         4  144° 


i6\        4          1440  4          24  4  V         4          24 


256      V       12         1440     }       10  V       12         1440     /\        4 
4 


1        f  \ 

Jo    1256  \ 

'    O" 

+   (  2  +  B'ijr 

4    \ 

7.  £ 
24 

f  (  &    ( 
=  Jol^6(10 

+-  -  — 
3 

/36   /s      p2 

£4\ 

6  MINIMUM  ERROR 

Putting  the  differential  coefficient  of  this  with  respect  to  b  equal  to  zero  we 
find  for  the  value  of  b  to  make  this  expression  a  minimum 

^_^+      0" 

8       i 1520 ' 

and  consequently  A  --  i  —       +       __  ^    , 

which  will  be  found  to  agree  with  the  expansion  of  2  cot2  £/3log  sec  A/3  as 
found  above. 

Substituting  this  value  of />  in  the  expression  for  M  we  find 

M*    I     /.••_____ 

768   1536   11520 

J?_   f^_    jS1" 
256   1536   11520" 

& 

,     to  3  terms  expanded, 
3  "4 

dividing  J/  by   i-cos/3  =      (i-     -  +  ...J   we  find  p.-  and  consequently 
a  -  -  (  i  +        ...  )  =        nearly.    If  we  substitute  the  value  found  for  b 

1/192  \       24/14 


in  the  original  expansions  we  find  for  the  scale  values  at  the  centre, 
where  9  -  o,  (dr;</0)0  -  (r/s'm  $)„  =  i  +  -- b  -  \  --  +  ...;  and  for  the  limiting 
radius  and  scale  values  thereat 

24      480      8960 

(r/sin  0)3  =  i  +  ^"+  ^V... 
8       96 

Supposing  instead  of  using  the  value  of  the  constant  A  which  makes 
the  mean  square  error  a  minimum,  we  had  used  A  ----  i,  as  did  Airy,  then 
we  shall  find  the  values  of  the  corresponding  quantities  by  making  b  -  o    » 
in  our  formulae  :  thus  (drfd6)n     (r/sin  6),,  -  r, 

ft'2        84^ 
48  '   320 


16       192 

...  )  :       ii          /  (        -  )  /S2  (  i  •)-  -     —  )  =        nearly  : 
i  o       /  \'    \3o4/  1 20 /        9 

so  that   both  the  maximum  and   mean  scale  errors  are  about  50  per  cent, 
higher.     Airy's  solution  with  A      \  does  not  even  give  the  minimum  error 
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solution  with  the  stipulation  that  the  scales  shall  be  true  at  the  centre,  as  will 
be  proved  later,  which  shews  the  importance  of  correctly  evaluating  the 
constant  A. 

To  find  for  what  value  of  0  in  the  Minimum  Error  solution  the  perpen- 
dicular scale  error  vanishes  put  r/sin  6  —  o  and  we  find  -  —  =  o  or 

1 6  1 6 

6  =  —r-  ft  nearly. 
V3 

It  might  appear  that  the  expansion  for  M  contains  no  terms  beyond  ft"; 
there  is,  however,  a  term  in  ft1'2  as  follows  : 

^_ 
30240 

Breusing's  Projection. 

As  we  have  already  remarked,  Breusing  conceived  a  projection  which 
is  the  geometrical  mean  between  the  Orthomorphic  and  Equal  Area,  viz., 
r—  2  v/(sin  i^tan-i^),  and  Hammer  has  pointed  out  how  nearly  this  agrees 
with  Airy's  solution  taking  A  \.  The  writer  finds  that  whether  the 
arithmetic,  geometric  or  harmonic  mean  of  these  two  functions  be  taken, 
we  obtain  a  projection  agreeing  with  Airy's  to  2  terms  in  the  expansion  of 
M,  and  to  at  least  that  order  in  the  other  quantities.  The  harmonic  mean 
leads  to  a  much  simpler  function  than  either  of  the  others  and  we  shall 
investigate  it  accordingly. 

The  Breusing  (Harmonic)  Minimum  Error  Projection. 

The  harmonic  mean  of  2  tan  \9  and  2  sin  \Q  is  4  tan  \0.  Consequently 
write  r^  ^A  tan  \6  and  determine  the  best  value  of  A  as  before.  We  shall 
find 

log  sec  -i-/3  -  tan-  -[ft  '  log  sec  \fi  —  tan-  [ft  ' 

To  shew  how  nearly  this  projection  agrees  with  Airy's  we  shall  find  in 
applying  it  to  the  extreme  case  taken  by  Clarke,  viz.,  ^=113°  30',  the 
limiting  radius  comes  out  as  1*5781  as  against  1-5760,  and  the  value  of  M 
as  0-1579  as  against  0-1569,  i.e.  its  average  scale  error  is  only  i  part  in 
3000  greater. 

By  expanding  it  in  the  same  way  as  Airy's  we  shall  find 


256  )       3      12      5760 


For  the  scale  values  at  the  centre  we  get 


At  the  limiting  radius  we  find 

B-       B4         6* 


/        B-       B4         6*        \ 
/?  (  i  -  —      —  --         -  .  .  .  J  ,  agreeing  to  4  terms  with  Airy's, 

R* 

- 


13-      /34 
(r/sin  0)0  =         i  +   8  +        ...  „          ,,3 

This  much  simpler  formula  thus  leads  to  almost  identically  the  same 
result  as  Aiiy's,  and  when  we  remember  that  in  the  extreme  case  of 
ft  1  13°  30',  a  map  covering  more  than  two-thirds  the  surface  of  the  globe, 
the  maximum  error  is  actually  72  °/o,  it  certainly  does  not  seem  worth  while 
using  the  more  complex  formula  of  Airy  to  reduce  the  mean  error  by  only 
i/3OOOth  part. 


The  Orthomorphic  Minimum  Error  Projection. 

In  this  \ve  have  r  =  2  A  tan  J,0  and  we  shall  find  that  A  has  exactly  the 
same  value  as  in  Airy's  and  J/  double  that  value,  viz., 


A  -  2  cot-  1(3  log  sec  i/3  =  i  -      + 


'-"""*   '"'  W 


11520 


J/=  4  sinj  1/3 

tan-/,/}  192  '  N/9&  \         24 

The  scale  value  at  the  centre  is 

(dr  </#)„  --  (r  sin  $)„  =  y/  =  i  -       + 

8       11520 

and  the  limiting  radius  and  scale  value  thereat  have  the  same  values  as  in 
Airy's. 

The  Equal  Area  Minimum  Error  Projection. 

In  this  we  have  r   -  2.-!  sin  .!#  and  on  intesratinti  we  shall  find 


/?< 


thf  term  ,.  being  the  diminution  in  .)/  caused  bv  the  use  of  the  con- 

^  ;o 

slant  A.    Tins  system   is  therefore  capable  of  practically  no  improvement 
bv  the  use  of  a  scale  fonstant. 
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The  scale  values  at  the  centre  are  A  =  i  -  -^     .,  and  at  the  outer  limit 

128 

of  the  map  we  have 


_  F  _  _£!_ 
8       192  •"' 

(r/sin*),=       !  +  £+-£.. 

The  use  of  a  scale  constant  with  2  sin  \6  is  in  fact  not  appropriate  to 
this  Projection,  for,  as  we  see,  it  destroys  the  characteristic  property  of 
reciprocity  of  scale  values,  at  any  rate  after  the  first  2  terms.  The  same 
remark  applies  to  its  prototype  r  ~-  4  cot  ^0  log  sec  \Q  for  which  also 

ft* 

M=         .    It  will  be  shewn  later  that  for  this  we  should  write  the  full 
192 

formula  r  ~  4  cot  \Q  log  sec  W  +  2B  cot  \A,  and  after  determining  B  so  as  to 
make  M  a  minimum  for  a  zone  from  0-a  to  @  —  (3,  we  put  a  =  o,  which 
amounts  to  making  B  =  Q  originally. 


The  Equidistant  Minimum  Error  Projection. 

Clarke  in  his  Paper  on  Airy's  Projection  published  in  the  Philosophical 
Magazine  for  April  1862  compares  it  with  two  others.  One,  the  Perspective 
Projection  of  James,  has  received  considerable  notice  but  the  other,  the 
Equidistant,  or  Equal  Radial  as  Clarke  calls  it,  is  seldom  mentioned, 
although  as  he  shews  when  drawn  with  the  same  limiting  radius  as  Airy's 
its  M  is  to  Airy's  as  1047  :  ioo'o.  Clarke  does  not  explain  very  clearly 
why  he  draws  the  limiting  radius  equal  to  Airy's,  but  what  he  does  has  the 
effect  practically  of  writing  the  equation  of  the  Equidistant  as  r~A9  and 
then  of  determining  A  so  that  the  mean  square  scale  error  shall  be  a 
minimum. 

On  applying  our  method  we  find 


A~ 


24      2160      36.',  880 


•  r  •     .>  i    rt  \  *•        *»»"•  r»    •    a  r»  /  Or*       I  P  4A-* 

^"-i"-,  _cos/3+  iVcosec^r  -^'  -»-  ^ 

Jo 

For  /?=  113°  30'  this  gives  ^=0-1634  as  against  0-1643  found  by  Clarke 
by  mechanical  integration,  for  the  case  r?  the  same  as  in  Airy's  solution. 
The  scale  values  at  the  centre  of  the  map  are 


io 


and  at  the  outer  boundary  we  find 

r&  =  B  ( i  —      -    -  ...},  agreeing  with  Airy's  to  2  terms, 

\  24         2IOO         / 

R-         7  /34 

(drldQ)*-        i  —        —--—...  ,,  „  i  term, 

24      2100 

(r/sin  0)0-          I  +    Q   4-  -!—-   ...  „  „  2  terms. 

8         i  oo 

The  value  of  the  mean  square  error  /i 

,    7/8'J 


This  exceeds  the  value  in  Airy's  in  the  ratio  N/J92  :  N/i72'8  =  105  :  100. 
In  the  writers  opinion  this  slight  increase  in  the  value  of  the  average  scale 
error  is  far  out-weighed  by  the  valuable  property  of  uniform  radial  scale 
and  the  simplicity  of  the  formula  of  the  Equidistant  Zenithal  Projection. 

If  we  make  A  =  i  we  shall  find  M  =  -~  as  against  -  -  in  Airy's 

216  768       153-6 

formula.  Therefore  Airy's  Projection  as  exhibited  by  the  author  and  other 
writers  is  not  per  se  so  good  as  the  Equidistant  in  the  ratio 

v/2i6  :  s/i53'6  =  1 19  :  100. 

The  merits  of  this  projection  in  comparison  with  others  are  described 
in  a  Paper  by  General  de  Coatpont  in  the  Bulletin  dc  la  Soriete  de  Geo- 
graphic for  1894,  but  without  the  use  of  the  scale  constant.  Its  invention 
is  attributed  to  Guillaume  Postel  in  1581. 

Modification  of  Minimum  Error  Projections  by  making 
total  area  of  map  true. 

It  will  be  observed  that  in  all  the  Minimum  Error  Projections  considered 
the  formulae  for  the  limiting  radius  r$  when  expanded  agree  to  the  first 

two  terms,  vi/.,  r$  --  ft  (  i  -  ...j.  Now  in  the  Equal  Area  Projection 
rp  ft  (  i  -  --...).  It  therefore  appears  that  we  shall  sacrifice 


very  little  in  the  value  of  the  average  scale  error  if  we  choose  our  scale 
constant  so  that  equality  of  total  area  is  absolutely  ensured,  and  this  is 
a  valuable  quality  in  any  map. 

In  some  cases,  as  we  shall  shew,  the  formula  for  the  quantity  A  is 
simplified.  In  every  case  considered  the  expansion  for  JM  is  altered  only 
in  the  term  (V"  and  upwards. 
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(i)   Airy's  Projection  with  total  area  true. 

Here  r$  =  A  tan  \ r  /3  +  2  cot  \  /3  log  sec  J  ft ; 

therefore  to  make  total  area  true  equate  this  to  2  sin  A/3,  whence 

A  =  2  cos  i/3-2  cot2  A /3  log  sec  l£=  x  -£+-£. 

8        192        7680 
Substituting  this  we  find 

._      (?  (         /34\  /83    /        /?-• 

J/  =  -V  I  i  +  ~-    ;    /*  -  -  -    i  +  -' 
384  V        384/  H92\        24 

= 2  sin  I A 


24         I920 

O/r/,/*),  =  i  +  £  ;     (''/sin  0),  =  i  +£!+  ||i  +...  =  sec  ^. 

3°4  o        384 

(ii)    Breusing  (Harmonic)  Projection  with  total  area  true. 
Here  r$  -  4  A  tan  \(3  -  2  sin  A/3  -  4  sin  ]  /S  cos  |/3, 

/3'2        yS4  5" 

whence  ^4  -  cos2  1/8  =  i  -    -  +  J  -  +  ..., 

16      768      92160 

£6  /        u/34\ 

and  J/=  -^-    i  +  -     -  )  j 

384  V        384o/ 

^  =  /3(i-^  +  -^...N)=:2sinA/3, 
V       24     1920     / 

(drldPfc  =  i  ;     (rlsin  6)ft  -  i  +  ?  +  ^  .  .  .  =  sec  A  /?. 

(iii)    Orthomorphic  Projection  with  total  area  true. 
Here  r$  =  2  A  tan  A  /3  =  2  sin  A  /3, 

)S-       /S4          /3fi 

therefore  A  =  cos  A/S  =  i  -  ~  +  —  ---  +  .  ..  , 

8       384      46080 


whence  ^ 

19 

-  (/"/'sin 


(iv)    The  Equal  Area  Projection  of  course  remains  unchanged  and 

we  have 

r$  —  2  sin  A/3,     A  —  i, 

i6  +  768J;    ^s'V)^1  ~  48'")' 
i 


(r/sin#)3=i+f'+^...  =  secA/3. 
o        .^04 


12 


MINIMUM  ERROR 


(v)    Equidistant  Projection  with  total  area  true. 
Here  r$  =  A/3  -  2  sin  i  /?, 

2    .  /?J        /34 

whence     A  =  -^sm-k8=i-    -+-         — ..., 
/j  24      1920 


„_ 


io,o8o,ooo/ 

fl"  /34  2 

(d>'l(16)&  =  i  —    -  H —         -...=:     sin  *B, 
24      1920  p 

(r/sin  0)^  =  i  +       +  ^—    +  ...  =  sec  i/?. 
The  perpendicular  scale  error  vanishes  when  0  -  ^fi  nearly. 

The  Minimum  Error  Perspective  Projection. 

The  geometrical  meaning  of  the  Perspective  Projection  is  as  follows  : 
Let  the  circle  in   Fig.    i   with  centre   O  represent  a  section  through   the 

P 

A' 

A 


centre  of  the  Earth,  and  draw  a  line  through  O  and  a  point  J/  on  the 
surface  which  corresponds  with  the  centre  of  the  map.  Then  choose  a 
point  /-*  on  the  diameter  MO,  prolonged  if  necessary,  from  which  points 
on  the  surface  of  the  Karth  can  he  projected  by  straight  lines  on  to  a  plane 
of  projection  //>  which  is  at  right  angles  to  the  line  MOP.  Let  A  be  such 
a  point  and  A'  its  projection,  then  the  distance  J\rA'  is  r  in  our  formula, 
and  the  angle  NO  A  is  ti  ;  denote  the  distance  OP  by  //.  and  the  distance 
NP  by  k.  The  radius  of  the  Earth  OA  is  taken  as  unity.  Draw  AQ 
perpendicular  to  MOly  and  let  N  be  the  point  where  the  line  MOP  cuts 
the  plane  />/>.  Then 


A'N 
AQ 


ArP 

Q.P 


and    A'N 


OP+  QO     /i  +  cosO' 
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which  is  the  characteristic  formula  of  this  class.  It  should  be  noted  that 
both  Germain  and  Thomas  Craig  make  the  mistake  of  saying  that  k  is  the 
distance  of  the  plane  of  projection  from  the  centre,  or  rather  Craig  says, 
"...making  h  the  distance  of  the  eye  from  the  centre  and  k  the  distance  of 
the  plane  of  projection  from  the  same  point...,"  whereas  k  is  the  distance 
of  the  plane  of  projection  from  the  eye  (i.e.  from  point  P). 

If  the  point  P  is  taken  at  P'  on  the  surface  of  the  Earth  and  the  plane 
pp  tangent  at  M  on  the  opposite  surface,  then  h  =  i  and  k  =  2,  so  that 

2  sin  0 
r-  =  2  tan  \  0, 

I  +  COS  V 

and  we  get  the  ordinary  Orthomorphic  or  Stereographic  Projection. 

If  we  take  h  -  o  and  k  =  i  we  get  r  =  tan  0  the  Gnomonic  Projection. 

If  we  take  k—  i  +  h  so  that  the  plane  of  projection  becomes  tangent  at 
M,  and  then  make  h  increase  without  limit  so  that  the  projecting  rays 
become  parallel,  we  get  r  -  sin  0,  the  Orthographic  Projection. 

Several  Perspective  Projections  have  been  put  forward  from  time  to 
time  with  the  plane  of  projection  tangent  at  M,  or  taking  k  —  i  +  h,  with 
values  of  h  ranging  from  i  to  2*73,  e.g.  those  of  Parent,  Lowry,  La  Hire,  etc. 

General  Sir  H.  James,  Director  General  of  the  Ordnance  Survey,  in 
seeking  for  the  best  values  of  h  and  k  for  a  map  covering  two-thirds  of  the 
sphere,  or  /3=ii3°3o',  either  by  intuition  or  luck  hit  upon  the  value 
h  =  i  '5  and  he  took  the  plane  pp  at  the  same  distance  from  the  centre  of 
the  Earth  as  the  plane  of  the  tropics,  that  is,  he  made 


--sn  2  3   30  =  0-399..., 
and  therefore          k-NO+  OP=  1-5  +0-399=  1-899.... 

Fiorini  shews  (pp.  213-4)  that  with  these  values  of  //  and  k  James'  Projec- 
tion gives 

/•£=  1-5812,     M=  0-17  20,     (dr^tt)^     0-6292,     (r/sin  (9)0=  17243. 

Clarke,  by  taking  k  so  that  ^  =  1-5760,  the  same  as  in  Airy's  Projection, 
finds  J/=  0-1718,  which  differs  very  little  from  the  value  with  James'  values 
of  //  and  k  ;  but  on  finding  by  integration  the  values  of  //  and  k  which 
make  Ma.  minimum  Clarke  obtains  //  =  1-36763,  k-  1-66261,  r$^  1-5737, 
M^  0-1626,  (dr\d6}^--  0-8053,  (r/sin  0)0  -17160. 

Clarke  on  integrating  the  value  of  M  brings  it  into  the  form 

J/=  4sirr  /3-  -j, 
where 


and  where  m  =  h+i,     n^/i  +  cos  /3. 


In  the  article  on  Mathematical  Geography  in  the  Encyclopaedia  Britannica 
(Ninth  Edition)  Clarke  puts  these  expressions  into  a  more  concise  form, 
but  the  writer  has  found  the  above  expressions  easier  to  compute  with  a 
table  of  Squares,  Cubes  and  Reciprocals  like  Barlow's.  The  differential 
coefficient  of  the  above  expression  with  respect  to  h  is  unmanageable,  and 
Clarke  then  proceeds  to  find  by  interpolation  the  value  of  h  which  makes 
B~IA  a  maximum.  The  writer  in  checking  the  calculated  values  of  log  jB'2/A 
in  Clarke's  example  discovered  a  slight  error  in  the  value  for  //  =  i*39 
which  should  be  0*420712  instead  of  0-420665.  This  correction  leads  to 
uniformity  in  the  second  differences  and  gives  //  —  1-368  for  the  value  of  // 
corresponding  with  the  maximum. 

Perspective  Projection  expanded. 

/£sin  0 

By  expanding  the  formula  r  -  from  the  outset  the  writer  has 

h  -r  cos  0 

succeeded  in  finding  expressions  from  which  to  compute  the  best  values  of 
//  and  /•  directly,  as  follows  : 

Write  k  -  (h  +  i)  (i  —  />),  so  that  b  represents  the  distance  from  plane 
tangent  at  point  M  to  actual  plane  of  projection  //.    Also  write  c  =  - — 

and  expand  cos  B.    Then  we  shall  find 

f        6-       6*    f      i\      &'    f        c       i 

r  =  ( i  +  In  sin  0  { i  +  --c  +  --  c(c--,-\  +  —  c(c- +  - 

2         4     V      67       8     V        39° 

&-  /       i 
i  +      [c- 

2  V      3 

whence 

dr'jd6      i  +        (c  — )  + 

2    V       3/        4 

/  C        3#V       A      5#V  •      c       r  M 
-  l>  \  i  +  -      K  -  -     +        (c-  -  ~+         I  , 

2  V       3/        4   V         2      3o/J 

3  9°/  " 

,  f        &"        0*    f       i\      1 

-0\l   +-f+        C(C -}...}. 

2          4     V       6/      J 
Forming  the  (juantity 

J/      I '  1  (dr  d&  -  i  Y  +  (r  sin  6  -  i )-'!  sin  6d6 

Jo 

in  the  usual  way  and  integrating,  we  finally  arrive  at 

f_l\  * 

3/ 

5V-1  r~"~   «  +  M4/    "J 
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4l4V     4/      3\  f  2) 

84  /  io7<r2      146*:      2t\      ] 

-    4or--    -  +  -  -)-"t 

8  \  3  15        3°/      J 


8 
/S-/         2\      F  fur      ic       31  \       1 

+  —      2f-       )  +  -   +  -'  ...V  . 

2  V        3/       3  V   4        4       I20/      J 

Differentiating  this  with  respect  to  />  we  find  the  value  of  b  to  make 
a  minimum  is 

4/         .  \3      4      247          \288       360      9607 '"J" 
Substituting  this  we  find 


M  =    - 


24  4  4 


.m/  , 

P  I 


5  24         40       240 


285<r3       !97^2      4°7^"       i 
24  40         480       20 

7£.a,3f 

4        8 


40  DO  400 

where  the  first  two  lines  represent  the  value  of  M  when  b-o,  or  the  pro- 
jection is  made  on  the  tangent  plane  as  supposed  by  Tissot  and  many 
writers  ;  the  second  two  lines  represent  the  decrease  in  M  due  to  moving 
the  plane  to  its  best  position  b  units  towards  the  centre  of  Earth  ;  and 
the  third  two  lines  represent  this  best  value  of  M.  We  can  now  find  the 
best  value  of  c  on  either  supposition,  by  differentiating  the  first  or  third  lines 
with  respect  to  c,  and  equating  the  result  to  zero. 

Taking  the  general  case  first  we  get  for  a  minimum 


to-S  +  ff      f  -     -  +        +  ff 

\  2      8/          \  3       40      60      480 

•solving  this  by  successive  approximations  we  find 


... 

8      256      61440  8\        32      23040 

Substituting  this  value  of  c  we  find 


. 

384 
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which  differs  from  the  same  quantity  in  Airy's  Projection  by  the  term  in 
j810  only.    We  shall  also  find 
(drjd6\  =  (r/s'm  6\  =  i  -  1> 

&       c.B*          /3s 

=  i  —  —  —    J  ....  agreeing  with  Airy's  to  2  terms. 

16      1536       24576 

_£!_  ^        £" 

'   24    48o    8960 


and  we  also  find 


2O         288OO 
-     (A+l)(l-£) 


3  V         32       46080 
We  shall  now  check  these  formulae  by  computing  the  values  of  h,  k 
and  M  for  the  cases  considered  by  Clarke  and  Fiorini,  viz.,  (3—  113°  30' 
and  (3  =  90".    We  find 

//  —  1*367     as  against  1*368  h—  1*471       as  against  1-470 

/t=i'66o  ,,          1-663  k  =-  2-040  ,,          2*038 

Af=  0*1613          ,,         0*1626  Jlf=  0*03951          ,,         0*03959 

The  above  expansions  for  h  and  k  are  therefore  sufficiently  exact  to 
enable  us  to  calculate  the  constants  of  the  Minimum  Error  Perspective 
Projection  for  maps  of  the  largest  extent,  and  to  calculate  them  directly 
instead  of  by  the  longer  and  more  troublesome  method  hitherto  employed. 

Perspective  Projection  on  tangent  plane. 

If  the  projection  is  made  on  the  tangent  plane  we  have 

24  i  4        °       4 

21 C  7 


<3* 

5          24         40      240^ "')  * 


whence  we  shall  find  for  the  minimum 
3 


i  - 
3 


/(•--/•  -r  i  -  i  -        -  ....    Also  (dr  dB)u  =  (r  sin  0),,  -   i  : 

3    V         200      / 

R-  ->{: 

i  -       ...  ;    (r/sin  0)fl  —  i  + 

20  2( 
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These  properties  differ  considerably  from  those  of  Clarke's  Minimum  Error 
Perspective  Projection.  The  value  of  M  is  however  smaller  than  for  any 
other  Zenithal  Projection  considered,  with  scale  values  true  at  centre.  With 

(3=  -  we  find  h-  2-04,  which  does  not  differ  much  from  the  value  2-105  °f 

the  Third  Projection  of  Parent. 

From  the  general  formulae  just  investigated  we  can  find  the  properties 
of  the  other  well-known  Perspective  Projections  on  both  the  tangent  and 
Minimum  Error  planes.  For  instance 

h  =  o,     or  c=  i  gives  the  Gnomonic,        r  --     (i  —  /;)  tan  0, 
h  =  oc  ,  or  c  =  o          ,,         Orthographic,  r  --=     ( i  —  /')  sin  6, 
A  -  i,     or  c  —  -i-         ,,         Stereographic,  r  —  2(1  —  i>)  tan  .',0, 
the  appropriate  value  of  b  being  found  by  substituting  the  proper  value  of 
c  in  the  above  expression  for />,  if  the  projection  is  desired  on  the  minimum 
error  plane.    It  may  be  noticed  that  if  we  take  c—  *,  we  obtain  a  Perspec- 
tive Projection  which  is  approximately  Equidistant,  for  it  agrees  with  the 
latter  to  terms  in  &-  and  B-  inclusive.    And  as  the  areal  scale  is  to  two 

Cfr 

terms  expanded  drjd&.rj^mB=  i -f    -(4^-1)...  we  can  get  a  Perspective 

Projection  which  is  nearly  Equal  Area  by  taking  c'=\,  or  // -•  3.  In  this 
case  b  =  o  to  a  first  approximation,  a  result  pointed  out  by  J.  I.  Craig,  p.  58. 
Though  it  has  all  these  inherent  properties  the  Perspective  Projection 
has,  in  the  writer's  opinion,  been  over-rated  in  the  text-books  written  for 
practical  cartographers.  General  Sir  H.  James  extols  it  above  Airy's  because 
it  is  a  real,  or  purely  geometrical,,  projection.  But  as  we  have  seen,  it  is  in 
its  general  form  awkward  to  compute  both  for  plotting  and  for  finding  the 
scale  errors,  and  it  differs  inappreciably  from  Minimum  Error  Projections 
of  simpler  form  such  as  the  Breusing  (Harmonic)  or  the  Equidistant. 

The  Minimum  Error  Zenithal  Projection  with  scale 
true  at  centre. 

The  fact  that  Airy's  solution  with  the  constant  A  =  i  does  not  give  the 
Minimum  Error  Zenithal  Projection,  subject  to  the  condition  that  the  scale 
shall  be  true  at  the  centre,  caused  the  writer  to  try  and  find  why  the 
ordinary  process  of  the  Calculus  of  Variations  failed  in  this  instance.  After 
solving  the  Perspective  Projection  with  this  stipulation  it  seemed  clear  that 

the  requisite  function  when  expanded  would  be  of  the  form  r  —  0—       6'' ..., 

0  C 

because  the  Perspective  Projection  agreed  so  closely  with  Airy's  in  the 
general  case.  Also  if  we  assume  r=0  +  afr'...,  where  a  is  a  constant  to  be 

determined  by  making  M  a  minimum,  we  find  Af—      f  ioa'--t-    -+       J, 
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tion.     On  substituting  this  solution   in  the   fundamental  formula   of  the 
Calculus  of  Variations 

rft 
8  U=  (SyJ/  in  our  case)  =  /    K&rdQ  +  H$  —  H^, 

Jo 
dP 
where        K=N     M  +  -<    #*  =  «>*  (/>-  ...)  j    #„  =  &-,,  (1\  ...), 

the  equation  S£7=  o  was  found  to  hold,  though  the  usual  separate  conditions 
A'=  o  and  H$  —  //„  =  o  did  not  hold. 

The  writer  was  fortunate  enough  to  interest  Professor  J.  W.  Nicholson, 
K.  R.S.,  in  this  question,  and  with  his  assistance  the  fundamental  equation 
of  the  Calculus  of  Variations  was,  at  any  rate  partially,  reconciled  with  the 
above  indirect  solution  in  the  following  way.  We  have 

M      I     {(t/rjdd-if  +  (r,  sin  6-  i)aj  sin  Odd, 
Jo 

and  in  the  usual  notation  of  that  Calculus 


o 

Now  assume  /  can  be  expanded  in  the  form  r  =  a6  +  l>63  +  c65..,,  where 
a,  l>,  c  ...  are  constants  to  be  determined  so  as  to  fulfil  the  conditions 
prescribed.  Then  since  drjdd  —  i  when  9  -  o  we  have  a—\  and 

r=e  +  !>&•  +  &'.... 
Now 

A'  -  r  l''jcr  =  2  (r/sin  6  —  i  ),    and    P  =  c  /-"/  rp     2  sin  ^  (dr'^iO  -  i  ), 
also 

/  =  dr',dO      i  +  tfP  +  5^4  ...,      r/sin  B  =  \ 
and  we  can  write  8/--  S/^6:;  --  o^6r'....  8/>  -  387^-  +  5o<:-6'4  —    Substituting  these 


3       360;  <SJ 
For  this  to  vanish  for  every  variation  of/'  we  must  have 


o, 
4  3 


and  if  we  cm  say  that  /'  must  be  independent  of  r,  then  c  -- o  and 
1    ..  Jl  J 

()O         320  60 
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and  we  shall  find 


6  \4<D      looo/      240  4° 

This  is  certainly  smaller  than  the  value  of  M  found  in  the  Perspective 

Projection   with   scales   true  at  centre,   viz.,   -^—  -  (i  +  —  —  ...  )     hut  the 

240  V         200       / 

evaluation  of  c  and  higher  coefficients  seems  to  leave  something  to  be 
desired  as  a  matter  of  principle.  So  far  as  the  cartographer  is  concerned 

the  formula  r--0  —  —  63  will  be  accurate  enough  for  all  practical  purposes, 

or  if  desired  the  next  term  in  the  formula  of  the  Perspective  Minimum 
Error  Projection  with  scale  true  at  centre  may  be  used,  vi/.., 

#:i      2i#!/3-      F 
r  =  Q~—  +  -  -  .... 

60       4000       150 

This  must  give  a  result  very  close  to  the  absolute  minimum  of  average 
error.  The  theoretical  difficulty  is  mentioned  in  case  it  may  attract  the 
notice  of  mathematicians  interested  in  the  Calculus  of  Variations,  which 
would  appear  to  require  development  in  those  cases  where  conditions  are 
imposed  on  the  derivatives  of  the  function  sought,  at  the  limits. 

Method  of  evaluating  the  average  error. 

There  is  another  theoretical  point  which  deserves  notice.  It  will  have 
been  remarked  that  Airy  and  Clarke  integrate  the  squares  of  the  scale 
errors  over  the  surface  of  the  EartJi  corresponding  to  the  map  concerned. 
It  appeared  to  the  writer  to  be  more  natural  to  find  the  average  error  by 
integrating  this  quantity  over  the  surface  of  the  map  and  then  to  divide  by 
the  area  of  the  map  before  evaluating  the  minimum.  It  appears  that  when 
once  the  form  of  the  projection  has  been  fixed  it  is  as  easy  to  use  one 
method  as  the  other.  For  example  taking  the  Equidistant  Projection 
r—AO  =  (i-jrl>)9,  the  scale  errors  will  be  expressed  as  before  but  the 

element  of  zone  on  the  map  will  be  expressed  by  2irrdr=2trr       dO,  so 

that  we  multiply  the  squares  of  scale  errors  by  rdrldQ  instead  of  sin  6.  In 
the  present  case  rdr'jdB  =  A'-Q.  Our  integral  M  is  however  divided  by 


so  that  the  A-  in  numerator  and  denominator  cancel  out  and  we  simply 
have  to  multiply  by  6  instead  of  sin  6  before  integrating  and  differentiating 
for  />.  By  using  the  expansion  method  we  obtain  in  fact 


40          I4000/          12     \  90 
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/3-  (        19/8-      43/34\ 

whence  b  =       -  (  i  +  +  -     -    , 

24  \         i  80       7560; 


and  J/  =  * 

40          I4000/  576    \  150          22O, 


40   14000 


32i/3v  -3^Yi  +  «£+  iifijTflvi 

140007   I728\    180   226,8oo7J 


1728  V         36°  " 

A'S- 

and  since  the  area  of  the  map  =  we  have 

2 


864;          V          boo 
as  against  fl  =  ^(-5-)  .  / 

found  by  integrating  over  the  Earth.  Integration  over  the  map  therefore 
gives  a  result  for  /x.  greater  in  the  proportion  of  i  +  /3"/6o  :  i  in  this  case. 

The  maximum  scale  at  the  boundary  is  however   i  +  —  +  — -  as  against 

8        004 

1  +    o  +       J  'n  tne  otner  method,  i.e.  slightly  smaller.    Application  of  this 
o        1 08 

method  to  the  other  fixed  functions  leads  to  similar  changes  in  the  third 
term  of  the  scale  values,  but  when  we  come  to  apply  it  to  the  problem  of 
finding  the  minimum  error  function  the  question  becomes  much  more 
difficult,  as  the  minimum  value  of  the  quotient  of  two  integrals  is  required. 
With  some  valuable  assistance  from  Prof.  A.  K.  H.  Love,  K.  R.S.,  the  writer 
was  able  to  put  this  problem  into  the  form  of  a  differential  equation,  much 
more  complex  than  Airy's,  and  an  equation  involving  values  at  the  limits. 
It  was  not  found  possible  to  solve  the  differential  equation  in  terms  of 
known  functions,  but  it  was  solved  in  series  for  3  terms  and  gave  a  result 
agreeing  with  Airy's  solution  to  2  terms.  The  question  of  finding  the 
average  error  by  integrating  over  the  map  is  therefore  not  of  practical 
importance  to  the  cartographer,  but  it  is  mentioned  here  as  a  point  of 
theoretical  interest,  and  also  to  reassure  those  who  may  be  inclined  to 
object  to  solutions  differing  from  the  absolute  minimum  error  solutions  in 
the  third  term  of  the  expansion  of  the  functions,  by  the  possibility  that  the 
validity  of  those  third  terms  may  itself  be  a  matter  of  argument. 

Summary  of  results. 

The  Minimum  Krror  X.enithal  Projections  have  been  thoroughly  examined 
by  deducing  and  comparing  the  formulae  for  their  mean  square,  or  average, 
scale  errors.  It  is  shewn  that  for  maps  of  the  largest  extent  the  Airy, 
Hreu^ing  and  ('larke  Perspective  Projections  are  almost  identical.  It  is 
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shewn  that  the  minimum  error  condition  differs  only  in  terms  of  a  negligible 
order  from  the  projection  obtained  by  making  the  total  area  of  the  map  true. 

It  is  also  shewn  that  the  Equidistant  Projection  with  total  area  true  has 
an  average  error  only  about  5  "/,  greater  than  the  average  error  of  Airy's, 
while  the  maximum  error  is  the  same' in  both.  This  small  defect  is  far 
outbalanced  by  the  simplicity  of  its  formula  and  this  Projection  is  therefore 
recommended  to  cartographers  as  the  best  xenithal  projection  in  all  cases, 
except  where  orthomorphism  or  equality  of  area  or  some  other  special 
property  is  required. 

The  cases  where  it  is  desirable  to  use  Zenithal  Projections  rather  than 
Conical  (of  which  Zenithal  form  a  special  case)  will  be  considered  later. 


CHAPTER  II 

THE  MINIMUM  ERROR  CONICAL  PROJECTIONS 

We  have  already  stated  that  the  Zenithal  Projections  are  only  a  par- 
ticular case  of  the  Conical  Class,  in  which  the  map  instead  of  being  developed 
on  a  plane  is  developed  on  the  surface  of  a  cone  having  as  axis  the  diameter 
of  the  sphere  which  passes  through  the  centre  of  the  zenithal  map.  This 
cone  is  then  supposed  to  be  laid  out  flat.  The  only  difference  is  that  the 
projection  of  a  small  circle  on  the  sphere,  perpendicular  to  the  axis,  instead 
of  being  a  complete  circle  of  circumference  2-rrr  is  a  portion  of  a  complete 
circle  2~nr,  where  n  is  the  sine  of  half  the  angle  at  the  vertex  of  the  cone 
and  is  called  the  constant  of  the  cone.  The  radial  scale  is  as  before  dr/d@, 
but  the  perpendicular  scale  is  now  #77 sin  6. 

The  expression  for  M  becomes 

M=  !   VdB  =  j  {(dr.dQ  -  i)-  +  (nr^rn  6  -  i )'-}  sin  6dO, 

J  a  J  a 

and  the  expression  for  the  corresponding  area  on  the  sphere  is  as  before 
(cos  a  —  cos  ft),  leaving  the  constant  multiplier  2-R-  in  each. 
In  the  usual  notation  of  the  Calculus  of  Variations 

jV^c^/dr=2«(Hr/sin0-i),     F=cVjcp  =  2  smO(dr,d6-  i), 

(IF,  dO  =2  cos9  (drdB  -  i)  +  2  sin  Od~r/d&-, 
so  that  the  differential  equation  N-dP^d6~Q  becomes 


The  complete  solution  of  this  equation  is 

r  =  A  (tan  1  0)"  +  (cot  \  0)"  i  (tan  '•  0)"  dO  +  B  (cot  i  6)", 

which  at  once  reduces  to  Airy's  solution  for  the  Zenithal  Class  when  we 
make  n  i,  for  I  tan  !  6//0  -  2  log  sec  \  0.  The  constants  A  and  B  are 

determined  from  the  conditions  that  P=o  both  when  O^a  and  6  ft, 
therefore  dr,dO--  \  when  0  -  a  and  0  =  ft,  or  the  radial  scale  is  true  at  both 
boundaries  of  the  map.  These  two  equations  lead  to  the  following  values 
of  the  constants  : 

(tan  ],ft)-"  -  (tan  A 


-  (tan  },ftf  I  "  (tan  \0)" 

B  •"— 

-"  -(tan  .U)-" 
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and  the  limiting  radii  of  the  map  will  be 

2  (tan  J/3)»  (    (tan  },&)" dO  2  (tan  ?,a)»  f   (tan  >0)V0 

^  =      (tan  i-0)2"  -(tan"$a)2ir~  '     ^  ~      (tan "pf"  -  (tan  .i'af™  ' 
so  that  the  depth  of  the  map  r$  —  ra 


_(tan  \Pfn  -  (tan  £a)2w         (tan  \fi)n  +  (tan  \a)n  ' 
The   only  trouble   so   far   is   that   the   integral    /  (tan  \6)ndO  cannot  be 

expressed  in  known  functions  except  when  //  i.  It  can  be  evaluated  for 
any  given  value  of  «  in  series  or  by  mechanical  integration.  The  calcu- 
lations may  appear  rather  formidable  at  first  sight  but  they  would  not  in 
any  given  case  occupy  an  expert  computer  more  than  a  day  or  two,  which 
is  a  small  matter  compared  with  the  total  time  occupied  in  the  compilation 
of  a  map.  We  shall  shew,  however,  that  by  means  of  expansions  and 
simplifications  even  this  labour  can  be  avoided,  at  the  expense  of  almost 
imperceptible  increments  of  the  average  scale  error. 

We  have  still  the  constant  «  at  our  disposal  and  in  order  to  find  the 
value  of  this  which  makes  M  a  minimum,  it  will  be  convenient  to  consider 
the  Orthomorphic  and  Equal  Area  parts  of  our  solution  separately. 

Minimum  Error  Conical  Orthomorphic  Projection. 

The  general  formula  for  this  is  r  —  2  A  (tan  |0)",  where 

,0 

A  =  I    (tan  |0)V0/{(tan  i/3)2"  -  (tan  .la)2"}. 
For  simplicity  we  shall  write  henceforth 

I      for     /    (tan  ?,  0)V0. 

Substituting  this  in  the  expression  for  J/  formed  in  the  usual  way,  we  shall 
find  after  some  transformations 

'0\-  /  "~l 

{(tan  1,B)-'1  -  (tan  'a)-"|    . 


/ 
We  can  also  expand  the  integral   I     in  terms  of  the  middle  value  of  the 

function  and  the  limits  by  the  following  general  formula 
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In  our  case 

AX)  =  (tan  ]X)"  ;  f"  (X)  -  '"Vn'v8  *"  (ta°  ^  ' 

Sill      X 

.  //4  -  6//3  cos  x  +  1>r  cos-  x  +  4"'  -  5"  cos  X  ~  *  CQS:'  X  /tan  i    y^ 

Thus 


84      1  6/v4 


36o'  sin4x 

We  also  have 

r? 

(tan  i/3/J"  -  (tan  'a)-"  =  2«  I    (tan  i6>)-"  cosec  Bdti, 

and  expanding  this  in  the  same  manner  we  get 

2ii  (tan  !x)~"  j        8-   4«-  -  6;;  cos  x  +  i  -  cos'-  x 
sin  x         1        6  sin"  x 

1  6//4  -  Sow"  cos  x  +  i  oo/r  cos'J  x  +  4o;r  —  7  o/v  cos  x  —  3o;/  cos3x 
+  i8cos'2x+  cos4x+  5 


120 


Dividing  out  we  find 


fY/!(tani/3)--(tan  'a)-; 

A   /   / 


f       8"   2«2  —  4«  cos  x  +  i  +  cos'-'x 
-i-z-  -    , 

A  (        6  sm-x 

48  w1*--  -96«3x+52//2cos2x—36«2+8o«cosx—  r6w  cos:ix 


_  _ 

6°  '  sin4 


Differentiating   the   last  with   respect  to  //  to   find   the   turning  value  we 
get  the  equation 

_  8-  4>t  -  4  cos  x 
6         sin" 


whence  ;/     cos  x  t.o  a  first  approximation,  and  inserting  this  value  in  the 
second  term  we  find 


When  we  insert  this  value  of//  in  the  expression  for  J/ we  find 

:      •  (  rt"  5*  8-          £4          )  40:'sil 

J/      40  sin  Y  •,  i  ...  —  i  +      + 

6       120  6       72       j  45 

and  dividing  by  the  area,  /j.~      ,-.  o4. 
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When  the  map  extends  to  the  centre  or  when  a  ^  o  we  proceed  as 
follows.    We  have 


Now 

where  /  =  tan  1/3  :  therefore 

a/s  V- 
(cot  W" 
.  i  / 


-cos/2-  211  (  f  Y(cot  1/3)4  . 

\  Jo    /  ) 

rp  f  fl+n  /:,  +  n  f,  +  n 

/    =.2  \  -   ---  +  -     __       \ 
Jo          [i+n      3  +  n      5  +  //        "  J  ' 


_ 

H 


.( i  +  n)-      ( i  +  «)  (3  +  „)         V  ( i  +  „)  (5  +  „)      (3  4.  ny)  •  • 

.  ,-2  \         n  2/-  /  2  I          \         1 

~"  4*   1 7      .      . ..  ~  /       — r^ ^v  +«(,  ,,  .  +  -,—  I  ...}-, 

+  n)  )      ) 


Differentiating  the  part  within  the  brackets  with  respect  to  ;/  to  find  the 
turning  value  we  get 

10  —  ztr      (9  —  /r)  (i  +  n}- \       "I 
i^~    ..\->  "•"         i  ~      \i         J  •  ••  r  ~~  ®> 


from  which  by  successive  approximation  we  find 
tz        2^  t~        /4 


«=  i--  --...        -  -        -  ...  =     --        _ 

2X2/36  29  8       144 

which  is  very  close  to 


as  usually  assumed. 

Inserting  the  value  of  n  found  in  the  expression  for  M  we  get 

rs\-j  i 

cot  .Jfl)2" 


.. .  ,,   and   tt  =  , 

48        /       12       768  v/384 

which  is  just  half  the  value  of  fj.  in  the  zenithal  system.  Thus  the  proper 
choice  of  the  constant  //  has  enabled  us  to  find  a  projection  which  has  only 
half  the  average  error  of  a  similar  projection  in  that  system.  Our  map  is 
however  now  fan-shaped  instead  of  circular. 

Minimum  Error  Conical  (approximately)  Equal 
Area  Projection. 

Taking  now  the  second  part  of  our  general  solution  we  have 
r  =  2  (cot  10)"   ((tan  £0)  V0  +  zB  (cot  A0)", 

and  on  forming  the  expression  for  J/and  differentiating  for  the  value  of  B 
to  make  this  a  minimum  we  shall  find  we  must  take 

B  =  ((tan  ia)-"  lf  -  (tan  i/3f   H  /{(tan  W  -  (tan  .la)2"} 
I  yo  -o  '/ 
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as  before  in  the  general  case,  and  on  substituting  this  we  shall  find  as  before 
that  M  reduces  to 

2  [cos  a  -  cos  /3  -  2n  (  I  )  /{(tan  i-/3)2w  -  (tan  Ja)2w}~]  . 

[_  V.'a  /    /  J 

The  value  of  ;/  will  be  the  same  as  in  the  Orthomorphic  system  whether 
the  map  extends  from  the  centre  or  over  a  zone  a  to  /3,  i.e. 

„-,-£+  /J4. 

8        144 

in  the  first  case  and  =    cos  x  (  '  +  —  •••  )  ni  the  second. 

\       3°      / 

And  when  we  form  the  expression  for  M  in  the  general  case 
r  =  A  (tan  £0)"  +  (cot  '  0)"      /  (tan  ?,0)V0  + 


it  will   be  found  that  there  are  no  products   of  the  several   terms,   and 
therefore  the  same  values  of  «  will  apply  ;   also  the  value  of  M  will  be 

B6 
half  the  value  just  found,  viz.,  M-  for  map  extending  to  centre  and 

T53" 

,  ..     28"  sin  Y  .. 
M  —  -  for  a  zone. 

45 

It  therefore  appears  that  the  Minimum  Error  Conical  Projection  is 
simply  the  arithmetical  mean  of  the  Minimum  Error  Orthomorphic  and 
(approximately)  Equal  Area  Projections  and  that  all  three  give  a  map  of 
the  same  total  size  and  limiting  radii,  the  intermediate  radii  in  the  two 
constituents  being  arranged  to  give  orthomorphism  or  equality  of  area  as 
the  case  may  be  and  the  absolute  minimum  error  radii  being  the  mean  of 
the  two. 


Numerical  example  of  Minimum  Error  Map  of  a  Zone. 

We  will  now  illustrate  our  results  by  computing  the  constants  and  scale 
errors  for  a  map  extending  from  a  =  30°  to  /5  =  90°  from  the  pole  of 
projection.  In  this  case  \  —  6°  and  cos  x  —  -1-  5  ar)d  wt-'  shall  adopt  this 

value   for    ;/    instead    of   the    more   correct   value    cos  Y  (  r  +        ...   )    for 

V         3°       / 

simplicity,  and  also  because  the  term  in  S'J  does  not  affect  the  first  term 
in  the  expansion  for  J/. 

It  will  first  be  necessary  to  evaluate  the  function   /(tan  \0}"  ii6  and  this 

we  shall  do,  so  as  to  give  seven  equidistant  values  between  30°  and  90°,  by 

/""" 
Weddle's    Rule.     This   Rule  for  evaluating  the  definite  integral    I      u,dx 

from  seven  equidistant  ordinates  //„,  //,,  i/.2,  ...  it,-,  is 

I  ".,    d.\  !//„     +      II.  ^     4      //,     -(-     //,;     +     l/.j     -r      5      (//|      +     tf..     -     //.-,)[, 
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where  //  is  the  value  of  the  interval  between  abscissae.     It  is  absolutely 
accurate  if  5th  differences  are  constant,  and  errs  in  excess  by  only  J^A8?/,, 


if  6th  differences  are  constant.     Applying  this  to  evaluate  /(tan  \fffidQ  for 
intervals  of  10°  we  tabulate  tan  £-0  for  intervals  of  o°  50'  as  under  : 


0 

¥ 

tan  £0          « 

=  (tan  J0)* 

0 

o°    o' 

o"   o' 

O'OOOO 

O'OOOO 

I 

i  °  40' 

0°  50' 

•0145 

'1204 

2 

3°  20' 

1°  40' 

•0291 

•1706 

3 

5     o 

2°  30' 

•0437 

'2090 

'2090 

4 

6°  40' 

3°  20' 

•0582 

•2412 

5 

8°  20' 

4°  10' 

•0728 

•2698 

6 

10°    o' 

5°    °' 

•0875 

•2958 

Sum 

•9166 

•5992 

5  x  '5992 

2-9960 

—  x  3-9126=  r  1  738.  The  interval  //  =  i  c  40'=    -         c.m. 

10  34377 


Integral^—  =0-0341. 


Applying  this  method  to  successive  intervals  of  10°  we  get  the  following 
values  of  the  integral. 

Interval  I  (tan  \d]^  dO  Total  from  oj 


|(tan  i^ 


0      -    10 

0-0341 

IO 

0-0341 

10°  —  20° 

•0630 

20° 

•097  1 

20°  —  30° 

•0821 

30° 

•1792 

30°—  40° 

'0980 

40° 

•2772 

40°—  50° 

•1123 

50° 

•3895 

50°  —  60° 

•1259 

60° 

•5154 

60°  —  70° 

•1393 

70° 

•6547 

70°  —  80° 

•1529 

80° 

•8076 

80°  —  90' 

•1671 

90° 

•9747 

90°  —  100' 

•1824 

100' 

1-1571 

I  OO°  —  I  IO' 

•1994 

I  IO° 

1-3565 

I  10°  —  I20J 

•2188 

1  20° 

r5753 

re 
There  is  one  case  when  the  definite  integral        (tan  W}tldO  can  be  expressed 

Jo 

in  known  functions  and  that  is  when  0  —  - - ,  when  it  can  be  expressed  by 

2 

the  \lr  functions  of  Gauss  as  -\\lfi     — -  )  —  «/r  {  -     -   )  -  .     In  this  case  since 

2  (.'   V    4    /         \    4    /J 

//  =  |  it  becomes  i  \\l/  (J-)  -  ^  (;^)J  -  0-9749  as  against  0-9747  in  our  Table. 
As  this  is  only  computed  to  4  figures  and  depends  on  .several  additions 
the  agreement  is  close  enough.  Using  these  results  we  shall  find  A  =  i  -0870, 
B  =  o'i  129,  and  we  can  tabulate  our  results  as  follows  : 
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$ 


n     -       O     O 


ON    —       10  co 

co     O       -     ri 
ro    >o       10  O 


t^    O 

C§       ON 


$ 

ON 
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The  first  two  columns  are  obvious. 

The  third  contains  the  values  of  I  tan  \BdB  just  found  and  underneath 

re 
is  written  the  sum        +  B,  as  this  is  to  be  multiplied  by  (cot  \OY  as  shewn 

Jo 

in  the  fourth  column. 

Columns  5,  6  and  7  are  obvious.     Column  8  is  obtained  by  taking  the 
differences  of  columns  4  and  5  and  multiplying  the  result  by  //r/sin^  since 

drjdO-  i  =     H~  [A  (tan  ^)"  -  (cot  '0)»  (  /"'  +  //)}  . 

sin  (/    i  \  JQ  j) 

These  differences  are  set  under  the  figures  in  column  4. 

Column  9  contains  the  squares  of  the  scale  errors,  those  of  the  radial 
scale  being  set  under  those  of  the  parallel  scale. 

Columns  10  and  n  are  obvious. 

Summing  column  1  1  by  Weddle's  Rule  we  find 

0-008967  x  10 

M-  =0-0011:6. 

57-3° 

8K  • 

sin  y 
1  he  formula  M  =      -'—-gives  0-00151,  the  difference  in  the  last  figure 

being  probably  due  to  taking  n  -  cos  x  instead  of  cos  x  (i  +  -    )  .    Dividing 


by  the  area  of  the  zone  we  find  /A  =  0-042.  Thus  we  have  a  projection  in 
which  the  maximum  scale  error  is  +  12-5  per  cent,  on  the  upper  parallel, 
+  8-7  per  cent,  on  the  lower  parallel  and  —  4-5  per  cent,  on  the  mid-parallel, 
while  the  mean  square  scale  error  is  4-2  per  cent.  The  standard  parallels 
will  be  near  40°  and  80". 

Map  of  Hemisphere  in  Minimum  Error  Conical  System. 

We  will  now  calculate  the  particulars  of  the  map  of  a  hemisphere. 
Since  a  =  o  we  put  B  o  and  taking  n  -—  cos  i/3  =  '07071  we  find  A  -  0-8369. 
We  calculate  the  following  Table  (p.  30)  in  which  the  interval  has  been 
taken  as  15°,  in  order  to  apply  Weddle's  Rule. 

We  have  thus  a  projection  of  the  hemisphere  in  which  the  maximum 
and  mean  square  scale  errors  are  just  about  half  those  in  Airy's  Zenithal 
system.  There  is  one  standard  parallel  at  about  65°  from  the  pole.  The 
map  will  be  fan-shaped  with  a  gap  -  2vr  (i  —  n}  104°  about.  This  shape  of 
map  might  not  be  so  convenient  for  some  purposes  as  a  circular  map  of  a 
hemisphere. 

Simplification  of  formulae  by  making  radial  scale  true. 

An  inspection  of  the  scale  errors  in  the  Tables  will  shew  that  those 
perpendicular  to  greatly  exceed  those  along  the  radii,  particularly  in  the 
case  of  the  zone  map,  where  the  radial  scale  is  true  at  each  boundary.  We 
shall  therefore  simplify  our  formulae  and  sacrifice  very  little  in  accuracy  it 


co    O        n     — 


-  2       Q  Q* 

p      o  o      o 


0     0 


O     0         O 

50      o 


O 
O 
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we  assume  that  the  radial  scale  is  true  always,  and  then  determine  our 
constants  so  as  to  make  the  perpendicular  scale  error  squared  a  minimum. 
We  write  therefore  r—  a  +  6,  where  a  is  a  constant,  then 

drjdO  =  i,     nr's'm  6  ••=  n  (a  +  0)  sin  8, 

rp  ( 

M-  I    (nr/sm  0  -  i )'-'  sin  6d&  =  -  cos  a  -  cos  ft  -  2//  (ft  -  a)  (a  +  £(/?+  a)) 

J a  (. 

.,/    ,        tan  }>R  i'&  OdB        i&9~d9\\ 

+  ir    a-  log  -  +  2a        -. !    I        —  H  . 

0  tan  ia  Ja  sm0      Ja  sin  8/j 

Unfortunately  the  integrals  lOcosecOdO  and  lO-cosccOdf)  cannot  be  ex- 
pressed in  known  functions.  They  could  if  necessary  be  evaluated  by  their 
expansions  or  by  mechanical  integration,  but  calling  them  for  a  moment 
S1  and  S-2  and  putting  dM[dn  =  o,  we  find  the  turning  value 

a  +  j-  (ft  -  a) 


tan  .',-£ 

-  +  2ao,  +  o.. 

tan  .id 


Putting  next  dM Ida  =  o,  we  find 

(3  -  a  - 


n  log  (         ^ 
0  \tan  W 

Substituting  this  in  the  expression  for  n,  we  find  after  some  reduction 
i  (/3  +  a)  (log  tan  i/8  -  log  tan  i-a)  -  ^, 

n  —  (p  —  a.)  ,  " "~'7^~7\ —  "T7T""n —        — i     \         c*~>         ' 

A2  (log  tan  !>p  -  log  tan  §a)  —  Of 
whence        ^W  =  cos  a  -  cos  ft  -  H  (ft  —  a)  (a  +  i  (ft  +  a)). 

If  the  map  extends  to  the  centre  then  a--o  and  a  -  o.  The  quantity 
log  tan  |/3 /log  tan  | a  becomes  infinite,  so  dividing  the  numerator  and 
denominator  of  the  right  hand  side  of  above  expression  for  //  by  it,  we  find 


therefore 

^  /3'; 


n—  i  — 


12          2l6o         90720 

14         4320          181,440 


864  V1     60 


which  is  2/5ths  the  value  of  Mm  the  Equidistant  Zenithal  system. 

If  the  map  extends  over  a  zone  from  0  a  to  0  =  (3,  the  expressions  for 
n,  a  and  M  can  be  simplified  by  expanding  as  follows,  for  which  method 
the  writer  is  indebted  to  Dr  T.  J.  I'A.  Bromwich,  F.R.S.  Putting 
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let  i/sin  (x-c)  =  A0  +  A1t  +  A2-  +  A,-  +  ..., 

so  that 

1  cos  x        A      i  +  cos2  x        A      cos  x  (5  +  cos'2; 

sinx'  sin~x'  sin3  x  sin4  x 

A*  =  sir?"  ^5  +  l8  C°S"  X  +  C°S4  X^' 
and  so  on  ;  further  write 

•+5        df.  .   /  .          .    8-  o4 


L  -  .  0  .2 

sill  (x-«)  V  6 


.       /• 
J 


_5   sm(x-e)        3 


.7—  ,     . 

sm(x-e)       3    V  10            56 

28Jj/ 
then  n  =  yv^  -,.-.,  •    Substituting  and  reducing  we  find 

/        8-      \  /        8a\ 

«  =  cos  x     i  +     -".)'     a  =  tan  X  (  l  ~  -  }  ~  X' 

\        3°      /  N        5 ' 


=  2  sin  x  I  6  -  --  +  -  -    ...     -  ,   , 
6       120  6 

which  values  for  ;/  and  M  agree  with  those  found  for  the  Minimum  Error 
solution  so  far  as  our  expansions  go.  We  are  therefore  justified  in  rejecting 
the  rather  complex  formulae  of  that  solution  and  adopting  the  uniform 
radial  scale. 

We  will  now  illustrate  our  results  by  calculating  the  particulars  of  the 
same  /onal  map  as  before,  extending  from  30°  to  90"  from  the  pole  of 
projection.  Putting  8  =  30°  =  0*5236  c.m.  and  x  -  60°,  we  find  a---  0-5899, 
n  --  0-5046,  and  our  Table  stands  as  under  : 


I 

2 

3 

4 

5 

6 

7 

8 

1)       ^ 

n  c.m. 
6  +  a 

sin  6 

wr/sintf 

w/sin  6  -  i 

Squares 
V 

/'sin 

e 

w/'sin  0  -  i 
in  Minimum 
Error 

System 

3"'       c 

5-5236 

0-5000 

1-1237 

+  •1237 

o'oi  5301 

0-0076 

50 

4-0-1252 

i 

•"35 

40 

"698  1 
•2880 

•6428 

r  *o  1  1  o 

+  *O  1  1  0 

121 

78 

+    '0130 

5°' 

•8727 
•4626 

•7660 

•9634 

-  -0366 

1340 

10 

26 

-     '0364 

60 

•0472 
•637  • 

•8660 

•9540 

—  '0460 

21  I  6 

18 

32 

-    -0452 

7° 

•2217 
•8  1  .  6 

"9397 

•9728 

-  -0272 

740 

6 

95 

—    '0280 

80 

•396  5 

•9848 

1-0178 

+  •0178 

317 

3 

12 

+      "0  I  60 

1  ;O 

'57o8 

1  'OOOO 

ro<)i2 

4-  '0902 

8136 

8  1 

36 

4-    -0870 
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The  columns  are  self  explanatory  and  in  the  last  are  given  the  perpen- 
dicular scale  errors  in  the  absolute  Minimum  Error  Projection,  for  com- 
parison with  those  in  column  5.  By  Weddle's  Rule  we  find  J/  =  0-00159 
as  against  0-00156  in  that  system.  The  formula  to  first  term  gives 
M=  0-0015  !•  We  see  that  the  map  has  two  standard  parallels  and  that  the 
scale  error  at  the  mid-parallel  60°  is 

/        8-\  /        8-\         8- 

•  (.  +  x)/»x-« -(.+-)(.--)  =  - 5 . 

To  find  approximately  where  the  standard  parallels  lie,  let  y  be  their  angle 
from  the  centre,  then  A;  (a  +  x  -  y)  =  sin  (x  -  y),  therefore 


cos  X     i  +  3Q      tan  X     i  -  ~     -  y    -  sin  X  (i  -  £/)  ~  7  cos  X, 

whence  y  =  +  —  8  =  ±  0-588, 

that  is,  the  standard  parallels  lie  at  approximately  one-fifth  the  depth  of  the 
map  from  the  upper  and  lower  edges. 

It  will  be  observed  that  the  mean  square  scale  error  is  independent  of 
the  distance  from  the  pole  of  projection,  or  what  we  may  call  the  mean 
co-latitude  of  the  map,  using  co-latitude  in  a  general  sense  as  angle  from 
the  pole.  This  is  an  important  result,  as  it  affords  us  great  liberty  of  choice 
in  fitting  any  particular  country  to  a  map  in  this  system.  We  have  only  to 
choose  our  pole  so  that  the  depth  of  the  map  28  shall  be  as  small  as 
possible,  and  as  the  choice  of  mean  co-latitude  x  enables  us  to  select 
almost  any  curvature  for  our  parallels,  we  can  suit  almost  any  shape  of 
country.  Hammer  has  published  Tables  to  facilitate  the  construction  of 
Conical  Projections  with  pole  in  any  position,  but  when  the  pole  has  a 
skew-position  it  will  be  found  that  generally  the  meridians  and  parallels  are 
lines  of  contrary  flexure  and  the  map  has  an  ugly  appearance.  Practically 
speaking  the  only  useful  positions  are  when  the  pole  corresponds  with  the 
terrestrial  poles  or  when  it  is  at  an  infinite  distance  in  the  perpendicular 
direction,  leading  to  transverse  cylindrical  projections. 


Minimum  Error  Plate-Carree  Projection. 

When  the  mid-parallel  in  the  Equidistant  Conical  Projection  is  of  co- 
latitude  90°  the  projection  becomes  cylindrical,  and  if  the  pole  of  projection 
is  at  infinity  in  the  direction  of  the  terrestrial  poles  it  is  called  the  Plate- 
Carree  Projection.  When,  however,  the  pole  is  at  infinity  in  the  transverse 
direction  the  projection  becomes  the  Cassini  Projection.  In  both  cases  can 
these  projections  be  improved  by  use  of  a  suitable  scale  constant. 

Considering  the  Plate-Carree  Projection,  let  us  now  count  the  angle  9 
from  the  equator  instead  of  from  the  pole.  The  quantity  n  becomes  xero. 
The  parallel  scale  error,  if  the  cylinder  is  tangent  at  the  equator,  is  repre- 
V.  i 
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sented  by  sec  0  —  i.    Let  now  A  be  a  scale  constant  and  let  us  determine 
A  so  that 

rS  V~ Az        i  +  sin  8  ~~1 

2   I    (A  seed  -  i)- cos  Odd  =  2\    —log  -„  —  2^48  +  sin  8 

./o  L  2          i  -  sin  8  J 

shall  be  a  minimum.    We  shall  find 

A  =  28/jlog  (i  +  sin  8)  -  log  (i  -sin  8)}  -  i  —  --  —        ...  : 

2  &;)  u 

M -  2  [sin  8  —  282/jlog  (i  +  sin  8)  —  log  (i  -  sin  8)]]  =  —      and    p.-     ,— . 

45  V45 

These  are  precisely  the  results  which  we  should  have  obtained  by  making 
X     90°  in  our  general  formula. 

If  we  take  A      i,  or  make  the  cylinder  tangent  at  the  equator, 

(  .  i  +  sin  6}        fr  8" 

'  i  •  i  »-\   A  i  A  _L    4_    I  /-\rr  _  •          ii 


-          -          J  log- 

&  i  -  sin  b) 


10 


l>y  creating  two  standard  parallels  instead  of  one,  we  have  introduced  a 

8- 

negative  scale  error  along  the  middle  parallel  of  -  .-,  we  have  reduced  the 

G 

mean  square  scale  error  in  the  ratio  ^45  :  ^/20  --  67  :  4-5,  and  the  maximum 
scale  at  the  extreme  parallels  has  been  reduced  from  sec  8  =  I  +—  ...  to 

8-\  8- 


i  -  z    =  i  + 

6, 


Minimum  Error  Cassini  Projection. 

Tissot  has  pointed  out  (p.  162)  that  the  Cassini  Projection  is  simply  the 
1'late-Carree  turned  through  a  right  an<zle,  so  that  what  was  the  equator 
becomes  the  principal  meridian.  Our  result  above  shews  that  by  reducing 
the  scale  along  this  meridian  by  one-third  the  scale  error  at  the  extreme 
east  or  west  edge  of  the  map,  creating  what  we  may  call  two  standard 
pseudo-meridians  (or  sub-parallels),  we  reduce  the  mean  square  error  of  the 
map  to  a  minimum.  A  device  of  this  sort  is  used  in  the  projection  of  the 
International  Map  on  the  scale  of  i/million.  The  scale  error  on  the 
principal  meridian  is  made  equal  and  opposite  to  that  on  the  extremes,  i.e. 

the  scale  constant  used  is  I  -    -  instead  of  i  -  —  • 

4  6 


Minimum  Error  Mercator  Projection. 

The  work  we  have  just  done  enables  us  to  find  the  Minimum  Error 
Mercator  Projection  which  is  simply  the  Cylindrical  Orthomorphic.    In  this 
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the  radial  or  meridian  scale  instead  of  being  uniform  is  the  same  as  the 
parallel  scale.    Using  (r—rv)  as  distance  from  equator,  we  have 


,  .  .      .  /7T         \  i  +  sin  0 

which  gives          (r  —  r*)  =  log  tan  I  -  +  \0  }  -~  !,  log  ---  =  —  -n  . 
•i  \4         /  i  —  sin  # 


The  expression  for  Afis  simply  double  that  in  the  Plate-Carree,  the  constant 
A  is  the  same,  but  must  be  applied  to  the  perpendicular  as  well  as  the 

x° 

.  8~  ~ 
45  474 


/ 

meridian  coordinates,  and  u  is  J2  times  as  great          / 

V 


Minimum  Error  Conformal  Projection  of  Gauss. 

As  the  Cassini  is  simply  the  Plate-Carre'e  turned  through  a  right  angle, 
so  the  Conformal  Projection  of  Gauss,  used  for  the  survey  of  Hanover  and 
more  recently  for  the  survey  of  Egypt,  is  simply  the  Mercator  applied  to  a 
meridian.  We  see  therefore  that  its  mean  square  scale  error  can  be  reduced, 

J-_>  J4 

by  multiplying  both  coordinates  by  a  factor  A  —  i  - ....  in  the  ratio 

6       72 

671  :  474  and  its  maximum  scale  error  in  the  ratio   3  :  2.    The  limiting 
depth  of  the  rectified  Mercator  or  width  of  the  rectified  Gauss  Conformal 

is  A  }>  log  ^  =  8,  the  same  as  in  the  Plate-Carree  or  Cassini-Soldner. 

0  i  -  sin  o 

We  shall  now  compare  our  Minimum  Error  Equidistant  Conical  Pro- 
jection with  some  others  of  the  same  class.  The  first  of  these  is  : 


Euler's  Map  for  Russia  with  two  standard  parallels. 

This  projection  is  described  in  detail  in  Germain,  pp.  179 — 183.  It  has 
the  property  that  the  absolute  errors  on  the  extreme  parallels  are  equal,  and 
also  equal  to  the  error  of  opposite  sign  on  the  mid-parallel.  This  leads  in 
our  notation  to  the  two  equations 

n  (a  +  x  ~  &)  ~~  sm  (x  "~  ^)  =  s'n  X  ~  n  (a  ~  X)  -  ll  (a  +  X  +  ^)  ~  s'n  (x  +  $)' 
whence 


The  map  extends  from  co-latitude  20°  to  co-latitude  50"  and  Euler  finds 
by  a  calculation  not  so  direct  as  ours  n  0-8098,  a  0^0853  4  53' 22". 
As  this  value  of  a  is  so  nearly  5'  -  0-0873  Euler  adopts  the  hitter  value  and 
so  fails  to  reproduce  one  of  his  conditions,  vi/.,  that  of  equality  ol  errors  at 
middle  and  extremes.  The  following  Tables  give  the  particulars  of  this 
map  calculated  from  the  Minimum  Error  formula,  from  Euler's  own  con- 
stants, and  from  the  correct  values  of  Euler's  constants. 
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Minimum  Error  Formulae. 
=  0-82156,     a  —  tan 


6      6  in  c.m.  and  r       sin  & 


20 

25° 
30° 
35° 
40° 

45" 
50° 


=  °'°7927- 


Mr/sin  9      wr/sin  6  -  i      F=  Squares  f  sin  6 

I'O286          +'0286          o'OOoSlS          0-000280 


4 
88 

134 


0-34907 
•42834 

•43633 
•51560 
•52360 
•60287 
•61086 
•69013 
•69813 
•77740 
•78540 
•86467 
•87266 

•95193 

A/"=  0-00003  20  by  Weddles  Rule 
=  o'oooo3i3  by  formula. 

Same  Map  with  Euler's  constants,  n  =  0-8098,  a  =  0-08727. 


44 
77 
26 

H 
332 


20 

0-34907 

0-34202 

1-0321 

+  •0321    0-001030    352  +  (o-oiio) 

•43634 

25° 

•43633 

•42262 

1-0032 

+  -0032 

10 

4 

•52360 

3°° 

•52360 

•50000 

•9895 

—  '0105 

I  10 

55 

•61087 

35'; 

•61086 

•57358 

•9857 

-•0143 

204 

1  17  -  (0-0082) 

•69813 

40° 

•69813 

•64279 

•9894 

—  '0  1  06 

I  12 

72 

•78540 

45" 

•78540 

•70711 

"9994 

—  -0006 

O 

o 

•87267 

5°J 

•87267 

•76604 

1-0147 

+  •0147 

216 

165  +(o'oi  12) 

•95993 

Jf=  0-0000352  by  Weddle's  Rule. 

If  Euler's  assumption  had  been  reproduced  we  should  have  the  scale  errors 
at  extreme  and  mid-parallels,  multiplied  by  respective  sines  of  0,  equal, 
whereas  we  have  the  figures  shewn  in  brackets  in  the  last  column. 

Same  Map  with  Euler's  correct  constants,  ;/  =  0*8098,  #  =  0*08535. 

0    0  in  c.m.  and  r     sin  6  nr^'mi)     ;/;-/sin  6  -  i      K=  Squares         Fsin  6 

20        0-34907        0-34202         1-0285         +0-0285        0'00o8l2        0-000278 +  (0-00976) 

•43442 

25  '43634 

•52169 
30          -52360 

35 
40 
45 
5° 
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The  assumptions  made  are  now  very  nearly  reproduced,  but,  as  we  see, 
Euler  by  using  erroneous  values  of  his  constants,  was  fortunate  enough  to 
reduce  his  mean  scale  error. 

To  calculate  the  formula  for  M  or  /A  we  should  strictly  speaking  substi- 
tute the  values  of  //  and  a  in  the  equation 

'     ~  -cosec 


M  =  cos  a  -  cos  /?  +  n1  (  ar  log  -  '     ~.     +  211  I    6  cosec  6dQ~-  I    fr-'cosec  QdQ\ 

0  tan  ha  Ja  Ja  ) 


- 


but  as  we  know  ft  is  to  the  first  term  independent  of  the  co-latitude  we  are 
entitled  to  take  the  value  for      -  -  and  write 


f  .    «  i  4-  sin  8)    . 

J/=  2  \  sin  o  -  2  Ad  +  },A-  loo;          .    ^  V  sin  Y 
1  to  i  -  sin  8J 


55  ,. 

=  20X1  —  —  +  -       -  ...  -  2A  +  hA'l  i  +  >-+-—'  sm  x- 
1         6       120  \         6       24/J 

Now  A  we  know  is  the  scale  at  mid-parallel 


Substituting  this  we  find 


120  V     \240 

Hinks,  pp.  83 — 87  (following  Close,  Textbook  of  Topographical  Survey- 
ing, p.  107  et  seq.),  considers  five  different  assumptions  for  this  projection 
including  Euler's  (Case  I).  We  will  now  exemplify  his  Case  III  in  which 
the  Scale  Errors  (not  the  absolute  errors)  along  extreme  parallels  are  made 
equal  and  opposite  to  that  along  the  mid-parallel.  The  maximum  negative 
scale  error  does  not  occur  exactly  on  the  mid-parallel,  except  when  \  =  goJ, 
but  it  is  very  near  to  it,  so  that  Hinks'  Case  V  is  practically  the  same  as 
his  Case  III.  We  have  then  the  following  two  equations 


sin(x-S)  sinx  sin(x  + 

whence  we    et 


a  +  x  ~  tan  x  •  8  cot  o        -  tan  x 

Scale  at  mid-parallel  =  2  tan  .18  cot  8  =  i  -  -  ...  the  same  as  in  Euler's 

4       24 

map ;  the  value  of  M  in  this  projection  will  therefore  be  the  same  as  in 

S ") 

Euler's  to  first  term  =  —     "  X  .    ^Ve  find  n     0-8238  :  a  =  0-0733,  and  Table 
1 20 

as  follows  : 
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e 
20° 

6  in  can.  and  ; 
0-3491 

-      sin^ 
0-3420 

nrjsin  0 
1-0174 

w/sintf-i       V—  Squares          Ksin  0 
+  0-0174         0-000303         0-000104 

•4224 

25° 

,    '4363 
•5096 

•4226 

"9934 

-    "0066 

43 

18 

3°° 

•5236 
•5969 

•5000 

•9834 

—    "0  1  66 

275 

137 

35° 

'6109 
•6842 

•5736 

•9826 

-    '0174 

303 

174 

40° 

•6981 

•6428 

•9886 

—    '0114 

130 

84 

7714 

45° 

7854 
•8587 

7071 

rooo3 

+    '0003 

0 

0 

5o° 

•8727 
•9460 

7660 

1-0173 

+    -0173 

299 

229 

J/=  0-000044  by  Weddle's  Rule 
=  0-000041  by  formula. 

The  position  of  the  standard  parallels  in  this,  and  in  all  cases  where 
the  scale  error  on  mid-parallel  is  equal  and  opposite  to  those  on  extreme 
parallels,  is  given  by  #/•-  \&,  i.e.  about  one-seventh  the  depth  of  the  map 
from  edges. 

We  will  now  discuss  two  remarkable  projections  of  this  class  invented 
by  the  Rev.  Patrick  Murdoch,  F.R.S.,  and  described  by  him  in  a  paper 
published  in  the  Philosophical  Transactions  of  the  Royal  Society  for  1758. 
But  first  we  will  .investigate  the  condition  that  the  total  area  of  the  map 
shall  be  equal  to  the  area  on  the  sphere  which  it  represents.  The  area  of  a 
zone  on  the  sphere  between  co-latitudes  a  and  (3  is 

2-rrR1  (cos  a  —  cos  /8)  =  2irR-  2  sin  8  sin  x 
The  corresponding  area  on  the  map  is 

mr  (r/s2  —  ra-}  =  ?nr  (r$  +  ra)  (r$  -•  ra)  =  2mrR- .  28  (a  +  x)- 
Therefore  for  equality  of  areas  /v8  (a  +  x>  -  sm  &  s'n  X>  or 

8"          84 

//  (a  +  v)/sin  Y  —  sin  8/8  =  i  —  -  +  -  —  — 

6  120 

Thus  the  Minimum  Error  Conical  Projection  has  to  a  first  approximation 
its  total  area  true,  as  in  the  zenithal  system. 

Murdoch's  First  Projection. 

Murdoch's  Paper  is  very  concise  and  he  does  not  give  the  full  proof  of 
his  result.  lie  takes  ;/  --=  cos  Y  simply,  and  then  gives  a  rule  for  drawing  the 
map  which  in  effect  amounts  to  taking  a  +  x  =  tan  x  sin  8/8,  and  this  we  see 
is  the  necessary  condition. 

Murdoch's  result  was  afterwards  taken  up  by  the  (ierman  mathematicians, 
Tobias  Mayer  and  Albers  (see  Germain,  p.  185).  The  former  professed  to 
have  discovered  some  Haw  in  Murdoch's  argument,  and  the  latter  elaborated 
the  result  and  also  developed  his  own  Equal  Area  Projection  of  this  class 
in  which  the  radial  scale  is  of  course  no  longer  uniform.  The  writer  has 
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examined  Murdoch's  original  Paper  and  can  discover  no  flaw  in  his  state- 
ments ;  and  it  is  remarkable  that  Murdoch  by  a  kind  of  intuition  hit  upon 
a  projection  which  has  to  a  first  approximation  the  same  mean  square  scale 
error  as  the  absolute  Minimum  Error  Projection  of  its  class. 


-  ~ 

For  if  we  substitute  A=  i--  +  --  in  the  formula 

6       120 


9  in  c.m.  and  r 
0-5236 
1-1304 

sin0 
0-5000 

wr/sin  9 
1-1304 

wr/sin  d  -  i 
+  0-1304 

F=  Squares 
0-017004 

Fsin  6 
0*008502 

•6981 

I  -3049 

•6428 

1-0150 

+    '0150 

225 

'45 

•8727 
i  "4795 

•7660 

•9658 

-    -0342 

I  170 

896 

1-0472 
1-6540 

•8660 

'955° 

•0450 

2025 

'754 

T22I7 
1-8285 

'9397 

•9728 

-    '0372 

1384 

1300 

1  1963 
2-0031 

•9848 

1-0170 

+    '0170 

289 

285 

I-5708 

roooo 

i  -0888 

+    -0888 

7885 

7885 

we  find  M  -        sin  Y  as  before. 

45 

To  exemplify  this  we  will  calculate  the  map  from  co-latitude  30''  to 
co-latitude  90°  on  Murdoch's  First  Projection.    We  shall  find 
n  =  cos  60°  =  0-5000,     a  —  tan  x  sin  &,&  —  x~  °'6o68. 


40 

50° 
60° 

70° 


90 

2-1776 

J/=  0-00164  by  Weddle's  Rule 
=  0-00151  by  formula. 

We  thus  find  scale  errors  practically  identical  with  those  of  the  absolute 
Minimum  Error  Projection  of  this  class,  and  a  mean  square  scale  error 
greater  in  the  ratio  ^164  :  ^156  only,  and  this  for  a  zone  map  of  a  very 
large  extent.  The  following  rather  disparaging  remark  of  Germain  (p.  186) 
is  therefore  not  justified :  "  Quelle  que  soit  la  methode  adoptee,  cette 
projection  ne  peut  convenir  pour  une  bande  de  plus  de  8  ou  10  degres  de 
latitude,  et  dans  des  limites  aussi  rapprochees  presque  toutes  les  projections 
coniques  ont  a  peu  pres  les  memes  avantages  ;  quant  a  la  question  de 
conserver  la  superficie  totale  de  la  carte,  elle  n'est  que  secondaire  et  oblige 
au  contraire  a  comprimer  les  superficies  et  les  distances  dans  la  /one 
moyenne,  tant  en  longueur  qu'en  largeur,  et  a  les  etendre  dans  les  deuv 
zones  extremes  dans  toutes  les  directions."  The  statement  as  to  the  com- 
pression and  extension  of  the  depth  of  the  map  is  certainly  not  correct,  as 
the  radial  scale  is  everywhere  uniform  and  true. 

Murdoch's  Paper  was  followed  a  month  later  by  a  letter  from  Mr  Wm. 
Mountaine,  F.R.S.,  on  "A  short  Dissertation  on  Maps  and  ('harts,"  in 
which  he  challenges  Murdoch  "to  draw  a  map  in  which  the  superficies  of 
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a  given  zone  shall  be  equal  to  the  zone  on  the  sphere,  while  at  the  same 
time  the  projection  from  the  centre  is  strictly  geometrical." 

Murdoch  replies  in  a  Note  printed  as  an  Addendum  to  Mountaine's 
Paper  giving  the  Rule  for  his  Second  Projection  which  is  a  Perspective 
Conical  Projection,  and  he  says.  "  But  here  the  degrees  of  latitude  towards 
the  middle  will  fall  short  of  their  just  quantity  and  at  the  extremities  exceed 
it :  which  hurts  the  eye."  He  then  indicates  the  Rule  for 

Murdoch's  Third  Projection. 

In  this  he  takes  a  +  x  ~-  tan  x  •  ^  c°t  S-  This  value  of  a  -  x  is  quoted  by 
both  Germain  and  Tissot,  but  it  does  not  appear  to  have  been  noticed  that 
the  condition  of  equality  of  total  area  exacts  a  corresponding  change  in  the 
value  of  //,  vi/., 


45 


Scale  at  mid-parallel  =r—  -  +  -  —  ...  •  M  --        sin  x  as  before.    Applying 

6       ] 20  45 

these  results  to  the  same  map  as  before  we  get  ;/  =  0-5265,  a  =  0-6451,  and 
our  Table  stands  as  under  : 


9      tfinc.rn.; 


40 

5°' 
6oc 
70° 

80° 
90° 


/ 


'754 
463 
6S6 

'0547 


U  e  thus  have  a  map  which  has  the  following  properties  : 

(i)     The  total  area  is  true. 

(ii)    The  radial  scale  is  everywhere  true. 

I'm)    Two  parallels  are  true  to  scale,  distant  about  1th  depth  from  edges. 

(iv)     I  he  scale  errors  at  extreme  parallels  are  equal  and  about  double 

the  scale  error  of  opposite  sign  at  mid-parallel. 
(\)     I  lit.-  average  scale  error  is  to  first  approximation  the  same  as  in  the 

Minimum  Knor  Projection  of  this  class. 
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(vi)    The  maximum  scale  error  is  less  by  'th  than  the  same  quantity  in 

the  latter. 
By  substituting  the  values  of  n  and  a  we  find 

o/o  "  2  O4 

Scale  at  extreme  parallels  =  tan  bo  =  i  +  —  +         .... 

3       15 

,,          mid-parallel 

This  Third  Projection  of  Murdoch  therefore  would  appear  to  be  on  the 
whole  the  very  best  of  the  Conical  Class,  better  for  all  practical  purposes 
than  the  absolute  Minimum  Error,  as  the  reduction  in  the  maximum  error 
more  than  compensates  for  the  slight  increase  in  average  error.  It  is  very 
easily  computed  whether  by  the  full  formulae  or  by  the  expansions  in  terms 
of  S-.  This  is  a  curious  reflection  on  Germain's  remark,  p.  189,  "  Ce 
troisieme  systeme  e'tant,  sous  tous  les  rapports,  inferieur  aux  deux  premiers, 
nous  ne  nous  y  arreterons  pas."  It  is  also  some  satisfaction  that  it  was 
invented  by  a  fellow-countryman,  who  though  born  and  educated  in  Scotland, 
was  Rector  of  Stradishall,  Suffolk,  England. 


Equidistant  Conical  Projection  with  total  area  true. 

As  in  the  Zenithal  System  the  condition  that  the  total  area  of  the  map 
shall  be  true  introduces  practically  no  difference.  It  affects  the  value  of 
M  only  in  the  third  term  of  its  expansion,  and  for  //  we  have  the  equation 


whence 


27r  I     n6dO  =  277(1  -cos/3), 
Jo 

»  =  2  (i  -  cos  /3)  /i-  -  4  sin-  J/S  /?J  -  i  -  ^"  -r   f 

12          l6o 


differing  in  the  third  term  only. 

Substituting  this  in  the  expression  for  J/  we  find 


864  V      60 
as  before. 

The  maximum  scale  error  at  the  boundary 

=  ?iB  sin  /?  -  2  tan  ';/3  8  =  i  +  -     +  -  —  ... 

12  120 

as  against 

7r\     :  ^  :  i3/^4_ 

36O/  12          2l()O 

The  following  Table  exhibits  the  map  of  the  hemisphere  in  this  system,  in 
which  ft  —  0797  : 
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e 

r 

sin0 

wr/sin  6 

wr/sin  0  -  i 

V—  Squares 

Fsin0 

0° 

O'OOOO 

O'OOOO 

0-7970 

-  -2030 

0-04121 

O'OOOOO 

'5° 

•2618 

•2588 

•8062 

-•1938 

•03756 

972 

3°' 

•5236 

•5000 

•8346 

-•1654 

•02736 

1368 

45° 

7854 

•7071 

•8852 

-•1148 

•01318 

932 

60° 

1-0472 

•8660 

•L/>39 

—  -036  1 

•00130 

113 

75° 

I  -3090 

•9659 

i  "080  1 

+  '  0801 

'00642 

620 

90^ 

1-5708 

I  "OOOO 

1-2519 

+  •2519 

•06345 

6345 

J/=o'oi68  by  Weddle's  Rule 
=  0-0167  by  formula. 

In  the  Minimum  Error  Projection  of  this  class  the  value  of  M  was 
0-0091  and  the  value  of  the  scale  error  at  /?  —  90°  was  +0-1835,  so  the 
merits  of  the  uniform'  radial  scale  are  not  so  near  those  of  the  absolute 
minimum,  in  so  large  a  map  as  a  hemisphere,  in  .this  class  as  in  the 
Zenithal.  It  is  still,  however,  much  superior  to  Airy's  for  the  hemisphere, 
in  which  J/was  0-0391  and  the  scale  error  at  the  boundary  was  0*3863. 

It  may  be  objected  to  this  fan-shaped  Conical  system  that  it  makes 
a  rather  ugly  map  when  n  is  greater  than  say  i,  as  there  is  a  wedge-shaped 
gap  =  (i  —  n)  360°  not  filled  up.  This  gap  may,  however,  be  filled  up  by 
continuing  the  map,  which  has  the  advantage  of  placing  in  juxtaposition 
those  parts  which  abut  on  the  boundaries,  so  that  continuous  measurements 
may  be  taken  across  those  regions.  Theoretically  speaking  this  advantage 
holds  however  small  the  limiting  radius  ft  of  the  map,  but  it  is  doubtful  if 
there  would  be  any  advantage  when  the  gap  is  less  than  30°,  or  ft  less  than 
about  60°,  for  all  bearings  from  the  centre  are  too  small  by  about  one- 
twelfth  and,  when  the  projection  is  so  near  a  complete  circle,  this  might 
be  inconvenient. 

Plate-Carree  and  Mercator  Projections  with  total  area  true. 

r/3 

In  these  cases  we  have  Aft  =        cos  Odd  =  sin  ft  and  therefore 

Jo 

A  =  sinpjp=i-^  +  -@- 

6        1 20 

as  against  i     fL  -  ." 

6       72 

in  the  absolute  minimum  system.      It  will  be  found  that  the  value  of  M  is 
not  altered  in  the  first  two  terms  of  its  expansion. 

Expansion  of  General  Formulae  in  case  of  Zone  Map. 

In  the  case  of  the  /.one  Map  the  limiting  radii  are  the  same,  whether 
the  map  be  Orthomorphic,  approximately  Equal  Area,  or  the  mean  of  the 
two,  vi/., 

2  (tan  .',/?)"   /'"  (tan  >0)"</fl  2  tan  .U  I     (tan  M)"dO 

',,  <  ~ 
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.u      f  f3  u        i  n\n  jn    f  (tan  ^)"  ~  (tan  -3>  a)" 

Depth  of  map  =  2  I     (tan  tj8)ndv  \ . 

/          \  £       I  l/f.,^>     1   Q  \'M          />„,-,     1    ..  \JJ1 

Now  we  have  found 

(tan  Ii6}nd0  =  28  (tan  |x)" 

S4    «4  —  6«3  cos  x  +  7«"  cos2  x  +  4^"  -  5^  cos  x  -  «  cos3 


-  +  .. 


r        84  /       ft- 

=  28  (tan  ^x)"  ] r T~  c°tJ  X  r  ^  7/     cos  x  (  i  + 

-28  (tan  ix)'1  1 1  —         cot'-'  xr  if  >l      cos  Y  simiAy. 
(.        120 


and 

We  have  also  found  that 


8'2  4«'2  -  6/1  cos  x  +  i  +  cos2 ; 
sin  x          I        6  sin2  x 


84  +  1  8  cos2  x  +  cos4  x  +  5 


120  sn 

and  replacing  2/1  by  n  we  find 


(tan  W  -  (tan  J«)»- 


sm 


;/4  —  i  on3  cos  x  +  i  °«2  •+•  2  5«2  cos'2  x  ~  35  w  cos  X  ~  x  5"  cos  '  X      1 
84  +  1  8  cos2  x  +  cos4  x  +  S 

+  ---- 

120  sirrx 

Dividing  out  we  find 
{(tan  ±PYn  -  (tan  ia)2»j/[(tan  Ift)"  -(tan  ia)"} 

_  (cot  ix)H  /       ?  3«2-3« 
2(6          sin2 

^4  —  9ow3  cos  x  - 


60 


,-> 

(cotM"  (T_84^_84c^x) 

-4 


i  —  — -  cot'2x  + 

60  24     J 


•\i-r~          ~  { 

(  4°      J 

and 


(cot  },v)n  f       84cot2x)    -f 

_-^'_  |i  +-         -H   if  «  -COSY  simpiy. 

2  {  24       j 


and  =  28  (  i  *  —  )  if  //  =  cos  x  Dimply. 

V  3°      / 

These  shew  how  nearly  the  depth  of  the  map  is  true. 
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We  shall  now  obtain  an  expression  for  the  radius  of  the  mid-parallel 
and  it  will  be  simpler  to  treat  the  orthomorphic  part  separately,  where 

r  =  2  A  (tan  \8f  =  |a  I    (tan  i0)"</0/{(tan  i/J)-'1*  -  (tan  }  a)-'1}     (tan  .V<9)" 

I—       J  a  J 

=  (using  above  expansions) 

8-  ^r  -  zn  cos  x  +  i  +  cos-  v]   , 
i  --  ^-  — ±\  (tan  16}" 

6  sin2  x  J 

=  tan  x  (  i )  (         j    )     if  H  --•  cos  x  (  i  +  -  —  )  • 

ri^i  /  O"' 

1  bus  when  v  =  v,  rv  =  tan  Y  I  I  - 

V       5 

the  same  as  in  the  Equidistant  System. 

Now   put   ^  =  x  +  €)  where    e    is    the  angle  from   mid-parallel    to   any 
other;  then 


=  (tan  !  v)»  ~  ~  cos'" 

"  '  sn:x 


=    (tan  .! -x)"  \  i  +  e  cot  x  +  £  cot  x  --  +  ^  cot  x  •  •  •  r   \t  n  --  cos  x  (  T  + 

I  3°      6  )  V       3° 

therefore 

S-j 


re  -  r.  --     i U     i  +        +    ....          e     i  --+--...,. 

\5/\3°6/         \        6      6      J 

This  gives  a  very  ready  method  of  calculating  and  plotting  the  Minimum 
Error  Orthomorphic  Projection  of  a  zone,  when  once  the  radius  of  the  mid- 
parallel  has  been  computed. 

In  the  approximate  Equal  Area  Projection 

r  =  2  (cot  I  6)n  i  (tan  1  <9)V<9  +  zB  cot  i-0, 


(cot  .',  a)'J"       (cot  v/3)"" 

we  shall  find  for  the  radius  of  mid-parallel 

r, 


L          (COt  \  a)'"     -  (cot  lfi)SH 

Expanding  this  expression  in  the  same  way  we  shall  find 

sin  x  /       8- 
r*        n     V      6 

and  evaluating  rc  —  rx  in  the  same  way  we  get  e  f  I  +-        ~  J  ;   so   that  the 

mean   of  this  and  the  orthomorphic  part  is  simply  -i\  —  rx      e,  the   Equi- 
distant value. 
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The  Albers  Equal  Area  Projection. 

The  Projection  we  have  just  considered  is  not  strictly  Equal  Area  since 
it  satisfies  the  differential  equation  drjdd-i  ^  —  (;/r/sin  0  -  i),  whereas  the 
strictly  Equal  Area  Projection  satisfies  the  differential  equation 

drjdO      sin  0/nr. 
The  complete  solution  of  the  latter  is 

:+-sitv.i#y 

H  / 

and  is  called  the  Albers  Projection.  When  the  Projection  extends  to  the 
centre  r=o  when  6  =  o,  and  r=  zsm^O/^Jn  simply.  We  shall  now  in- 
vestigate what  value  of  ;/  makes  the  mean  square  scale  error  a  minimum. 

Putting  >Jn  =  I  +  a,  so  that  —.--  -  i  —a  +  rr—  ...,  we  have 

r  =  (i -«  +  o^..)  («--+  —  ...), 

V       24      1920      / 


Q-       ~  64 

nr/s'm  8  =  i  +  -$  +  -  -+...+ 
8      384 

M=  f    {(drjdO  -i)'  +  (r/sin  0  -  i)-}  sin 
Jo 


JO      1L 

+  64V  '  487    '     4    V*  '  48 
whence  for  a  minimum  we  shall  find 


,- 

16      512 

^ 
or  I  +  a=I-!6  + 

/?• 
and  (i  +af  =  n=  i  - 


which  value  of  n  differs  very  little  from  cos  .!/J=  i  -  -Q     ...  proposed 

o       3  *r 

by  M.  Tissot  to  make  the  maximum  scale  error  as  small  as  possible.     The 

(B~         /J4        x 
l+L  and   the   maximum    scale   error 

48      7660      / 


10      1^30 
third  terms  found  for  the  approximate  Equal  Area  Minimum  Error  system. 

06 

The  value  for    I/ is    --  to  the  first  term,  the  same  as  in  the  latter  system. 
768 
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As  to   the   Albers   Projection   for  a  zone,   r  =  (  r*  +  -  sin2  ^0  J   ,    the 

writer  has  not  succeeded  in  deriving  by  rigorous  mathematical  processes 
the  values  of  the  constants  r0  and  n  which  render  its  mean  square  scale 
error  a  minimum,  the  requisite  integrations  and  expansions  being  very 
complex.  By  writing  the  formula  in  the  form  r  ---  {a  (c  -  cos  0)p,  where 
a  =  2/n  and  a  (c—  i)  =  r,,'J,  we  have 

drjdO  --  a  sin  0/2  {a  (c  —  cos  0)j  2,     nr/sm  &  ~  2  [a  (c  —  cos  @)}-/a  sin  8, 
and 


Jasm-6       4(r-cos(9)      4(V-cos<9)^ 
M=  4-  .  —  —    -  -f  2  sin 


We  can  now  expand  the  integrals  in  terms  of  the  middle  value 
X  =  I  (ft  +  a)  ar>d  the  limits  8  =  -1-  (/8  —  a),  and  we  may  put  in  the  expansions 

a=  2/H  --  2/cosx  (i  +<?),     <r-cosx  =  sin'-'x  (i  +/)  2  cos  x, 
where  ^  and  _/  are  corrections  depending  on  8-  and  higher  terms,  the  values 
of  which  are  to  be  found  so  as  to  make  M  a  minimum.     The  writer  has 
not   taken   the  expansions   beyond   8-  and   so   far   he   finds   the   equation 
<'+_/+  .\8-  =  o  must  be  satisfied. 

By  analogy  with  the  approximate  Equal  Area  function,  we  may  however 
assume  that  if  we  make  the  scales  at  upper  and  lower  limits  equal  and  the 
total  depth  of  the  map  true,  we  shall  not  be  far  from  the  absolute  minimum 
mean  square  error  conditions.  These  assumptions  have  been  worked  out 
by  Fiorini  (p.  269)  and  amount  to  taking 

i  +  cos  a  cos  (3,  2  8-  cot  8 


cos  x  cos  8  cos  x  sin 


C  = 


whence  C  —  cos  x  z • —•       ci  =  — : 

2  cos  x  cos  x 

thus  e  -+   ,  ,  /  = and  our  equation  <-+/+  .!  82  -  o  is  satisfied,     ^^*e  also 

6  2 

find 

rv  =  tan 


the  same  values  as  in  Murdoch's  Third  Projection. 

To  exemplify  this  result  we  shall  compute  the  same  map  from  a  =  30° 
to  (3  --  90°  (Upper  table).  We  find  n  -  0-5265,  a  =  37988,  C=  1-1547. 

Tissot  (pj).  222 — 6)  has  investigated  the  values  of  the  constants  of  the 
Albers  Projection  to  make  the  maximum  scale  error  as  small  as  possible. 
This  involves  making  the  scales  at  upper  and  lower  parallels  equal  and  also 
equal  in  value  to  the  maximum  scale  of  opposite  sign  which  occurs  near 
the  mid-parallel.  We  shall  now  compute  the  same  map  on  Tissot's  assump- 
tion, for  comparison  (Lower  table).  In  his  notation  r^  (;-„-  +  -sin-  \  0V" 

\          n          '    /  ' 
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cos  #  Squares 

6       C-COS0  a(C-COS0)     ;-         «/-/sin0       dr\dB      «>-/sintf-i    dr\dd  -  \    and   f       Fsinfl 

30°     0-8660       1-097       1-048       1-104       0-906       +'104        -'094       -oiio       -0099 
•2887  88 

•0198 

40°       -7660      1-477       1-215        '995        ro°5        -'005        +'005       -0000       -oooo 
•3887  o 

"OOOO 

50"       -6428       1-945       1-395        '959       i '043        -'°4i        +'043       'oo  1 6       -0026 
'Sll9  1 8 

•0034 

60°       -5000      2-487       1-577        -959       1-043       -'041        +'043       -0016       -0029 
•654?  18 

•0034 

70        -3420      3-087       1-757        -984       i -oi  6       --016       +-016       -0003       -0006 
•8127  3 

•0006 

80        "1736      3727       i'93i       1*032         -969       +'032        -'031       -ooio       -0020 
"9811  10 

"OO2O 
90'    "OOOO    4'386    2*094    I'IO2      "907    +'IO2     —'O93    "0104    "0190 

1-1547  86 

"0190 

By  \Veddle:s  Rule  J/=o'oo3ii  as  against  0-00300  for  double  Murdoch's 
Third  Projection  and  jn  =  o'o6o,  so  we  may  be  sure  we  are  very  close  to  the 
conditions  for  absolute  minimum  of  M. 


70 


90- 


1-1550 

—  'O/O 

49° 

0-01052 

•2340 

•9360  1-247 

•970 

I  '03  1    —  -030 

90 

1*555° 

+  '031 

c6 

786 

'357i 

1-4284   1-431 

'934 

I  -07  1   -  -066 

436 

2-0474 

+  •071 

504 

940 

•5000 

2'oooo   i  "6  1  9 

'935 

[  -069   -  -065 

422 

2-6190 

-t-  -069 

476 

898 

•6580 

2-6320   i  '803 

'959 

1-043   -"°4i 

1  68 

3-2510 

+  '043 

185 

353 

•8264 

3-3056   1-981 

roo6 

•994   +  -006 

4 

3-9246 

—  -006 

4 

8 

•oooo 

4"oooo   2"  149 

1-075 

•930   +'07  5 

562 

4-6190 

-  -070 

4?o 

1052 
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H=  maximum  positive  scale  value  near  mid-parallel,  8    =  co-latitude  of  that 
parallel  of  maximum  scale  error,  then  if 

8  =  ?,  ((3  —  a),    H=  (sec  8)^  =  1*075,     n~cosx^=^'     cos  8^  =  nHz  =  0-5773, 
8^-54°  44';     rt)  --  (i  -cos  8fjC)/nff-  07868,     ?y  =  0-6190. 

The  maximum  scale  error  has  been  reduced  from  10-3  per  cent,  to 
7-5  per  cent,  at  the  expense  of  increasing  the  mean  square  scale  error  from 
6-0  per  cent,  to  6-9  per  cent.  It  will  be  observed  however  that  this 
maximum  error  occurs  at  the  very  edges  of  the  map  and  then  rapidly 
diminishes,  whereas  in  Tissot's  system  the  errors  between  the  standard 
parallels  are  much  larger,  being  50  °/o  greater  at  the  mid-parallel. 

When  the  mid-parallel  in  Albers'  Projection  becomes  the  Equator  the 

formula    changes    into    r  —  rtl  —  —,  sin  (/>,   where  </>  =  90°  —  0  or  angle  from 
equator.     The  minimum  value  of  J/will  be  found  by  taking 

;/  =  sin  8/8  —  i  —  -    +  ... . 
6 


Map  of  the  World  in  three  parts  in  Minimum  Error  Conical 
System  as  proposed  by  M.  Collignon. 

According  to  Germain  (p.  103)  M.  Collignon  proposed  to  construct  a 
map  of  the  World  in  three  parts  by  projecting  the  portion  of  the  surface 
from  the  poles  to  latitude  30°  on  the  Isospheric  Stenoteric  System  of 
Lambert,  i.e.  the  Conical  Equal  Area  System  or  r=  2  sin  \B\^jn,  and  the 
middle  zone  covering  30°  on  each  side  of  the  Equator  on  the  Isocylindric 
System.  The  greatest  error  of  distance  is  said  to  be  15  per  cent,  and  of 
angle  8°  12'. 

If  we  use  the  Minimum  Error  Equidistant  Conical  System  instead,  we 
have  for  the  parallel  scale  error  at  60°  from  the  pole 

+  B-/I  2  —  (c.m.  of  6o°)2/i  2  =  9  per  cent., 

and  for  the  middle  part  we  have  in  Murdoch's  Third  Projection  parallel 
scale  error  at  30°  from  Equator-  +  o'2/3  =  (c.m.  of  30°)-,  3,  i.e.  the  same,  so 
that  the  sections  would  have  a  rolling  fit  along  the  boundaries.  The 

greatest  alteration  of  angle  would  be   2  sin-1  -  about    ;°.     The  mean 

2-09 

8- 
square  error  in  the  conical  parts  =  i  —       —0-91. 

This  is  probably  the  best  map  of  the  World  that  can  be  obtained  in 
three  separate  sheets.  The  middle  portion  of  the  map  covers  exactly  half 
the  total  surface  of  the  sphere  and  the  end  portions  each  one  quarter.  The 
scale  along  the  Equator  is  of  course  too  small  by  half  the  excess  at  latitude 
30°,  or  1  >y  about  4!  per  cent. 
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G.  W.  Hill's  Projection. 

The  celebrated  American  astronomer  and  mathematician,  (1.  \V.  Hill, 
in  a  Paper  published  in  the  Annals  of  At 'at 'hematics,  Second  Series,  Vol.  10, 
No.  i,  October  1908,  investigates  a  projection  which  he  claims  to  be  the 
best  possible  on  Tchebycheff's  (or  Chebychev's)  Principle  that  "A  map  is 
the  best  when,  within  its  limits,  the  logarithm  of  the  scale  has  the  least 
possible  oscillation."  This  Paper  is  rather  difficult,  but  it  appears  to  be 
stipulated  that  the  map  shall  be  conformal,  i.e.  orthomorphic,  and  the 
principle  leads  to  the  result  that  the  scale  error  must  vanish  all  round  the 
boundary. 

The  scale  along  the  parallels  is  no  longer  uniform  and  the  meridians  in 
addition  to  having  a  variable  scale  are  no  longer  straight  lines 

Hill  applies  his  results  to  an  example  extending  from  latitude  40"  to 
latitude  50°  and  for  a  difference  of  longitude  extending  5°  on  each  side  of 
a  central  meridian.  He  gives  a  table  of  scale  errors  of  which  the  sum  is 
0^0754  for  66  values,  the  arithmetic  mean  of  which  is  o-ooii4  and  the 
maximum  o-oo254. 

If  Airy's  Zenithal  System  had  been  applied  to  the  same  map  it  would 
be  comprised  in  a  radius  (3  of  about  6° -3,  and  so  the  mean  square  scale- 
error  fj.  would  be  /3-/x/i92—  0*000906,  the  average  error  =  o'yyyg/x  0-00072, 
and  the  maximum  error  (3-,  8  =  0-0018. 

If  the  Conical  System  were  applied  to  the  same  map  a  /.one  to"  wide 
wrould  cover  it,  the  mean  square  scale  error  would  beS2/6'7  (where  8-  c.m. 
of  5°)  =  o-ooi  i,  and  the  maximum  error  would  be  <r,  3  -  0-0025. 

It  does  not  appear  that  Hill's  Projection  has  any  advantage  over  those 
we  have  investigated  and  it  is  certainly  more  difficult  to  compute. 

It  still  remains  an  open  question,  perhaps,  as  to  whether  there  is  a  class 
of  projections,  i.e.  of  developable  surfaces,  better  than  the  Conical.  The 
writer's  impression  is  that  there  is  not,  but  so  far  he  has  not  been  able  to 
prove  this  absolutely. 

Best  projection  for  a  given  country. 

As  we  have  already  stated,  the  best  projection  for  a  given  country 
depends  on  the  shape  of  the  country.  Tissot  lias  developed  a  method  of 
finding  such  a  projection  on  the  principle  of  making  the  maximum  scale 
error  as  small  as  possible.  He  applies  this  to  finding  the  best  projection 
for  a  map  of  Spain  and  (p.  69)  he  gives  the  resulting  scale  errors.  The 
map  extends  from  latitude  36"  o'  to  latitude  43  50'.  The  representation 
appears  to  be  nearly  conformal.  The  maximum  scale  error  is  -;  30119 
and  the  mean  square  scale  error  works  out  at  C-QOI  16. 

If  the  Zenithal  System  had  been  used  it  appears  that  Spain  could  be 
included  in  a  radius  (3  of  about  5'  =  0-087  c.rn.    If  the  Kquidistant  Minimum 
Error  Projection  were  used  the  maximum  scale  error  would  be  (3-  t 
and  the  mean  square  scale  error  would  be  /iJT3'T     0-00058.    This  is  there  - 
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fore  better  than  Tissot's,  for  we  could  have  conformality  with  the  same 
maximum  error  and  a  mean  square  error  of  /3'-'/x/96  =  0-00077. 

If  however  Murdoch-'s  Third  Projection  were  used  we  could  take 
/3  —  a  =  8°  and  8  —  4°  —  0-070  c.m.  Hence  the  maximum  scale  error  would 
be  S-/3  =  °'OOI6  and  the  mean  square  scale  error  would  be  8-7^/45  ^0-00073. 
This  could  be  made  conformal  like  Tissot's  with  the  same  maximum  error 
and  by  increasing  the  mean  square  error  to  o  00073  x  \/2  =  0-00102. 

Either  of  these  projections  is  therefore  at  least  as  good  as  Tissot's  and 
both  are  easier  to  work  out. 

Best  projection  for  Italy. 

Mr  J.  I.  Craig,  of  the  Egyptian  Survey  Department,  in  his  pamphlet  on 
Map  Projections,  raises  the  question  of  the  best  projection  for  Italy,  a  long 
narrow  country  with  a  general  bearing  of  about  45°.  It  would  appear  that 
a  zone  about  5.!"  wide  would  cover  Italy  proper,  excluding  Sicily.  If  we 
used  an  oblique  conical  system  we  could  represent  Italy  on  Murdoch's 
Third  Projection  with  a  maximum  scale  error  of  (o'048)"/3  =  0-00077  and  a 
mean  square  scale  error  of  0-00230/13-1  =  o'oooi  7.  To  represent  the  same 
country  on  the  Equidistant  Zenithal  System  we  should  have  to  take  /3-about 
51,'"  and  then  the  maximum  error  would  be  (3r/S  =  o'ooi  15  and  the  mean 
square  error  would  be  fl-li^'i  -  0-00070. 

It  appears  therefore  that  .we  shall  only  do  better  in  the  Zenithal  System, 
for  countries  of  moderate  extent,  if  ($-]  13-1  is  smaller  than  82/6"j,  that  is  if 
/?  :  8  is  less  than  1-4  :  i,  or  if  the  country  to  be  represented  can  be  included 
in  a  parallelogram  whose  longer  to  shorter  side  is  less  than  about  i '4  to  i. 
It  must  be  remembered  that  the  Conical  is  always,  theoretically  speaking, 
better  than  the  Zenithal  System,  but  we  are  excluding  the  conical  projection 
up  to  pole  of  projection,  where  such  pole  would  occur  in  the  middle  of  a 
country  of  moderate  size  and  produce  a  fan-shaped  map  with  a  small  gap 
=  360°  (i  -//)=  15°^. 

As  we  have  remarked  Hammer  has  published  Tables  for  Conical  Pro- 
jections with  pole  in  any  position.  The  formulae  are  capable  of  a  little 
simplification  when  the  axis  is  inclined  at  45°  to  the  meridian  as  in  Italy, 
but  even  then  this  skew  system  cannot  be  recommended  to  cartographers, 
as  the  decrease  of  scale  errors  will  be  more  than  counter-balanced  by  the 
disadvantages  of  complex  formulae  and  contrary  flexure  of  the  projected 
meridians  and  parallels. 

Summary  of  Results. 

The  Minimum  Error  Projections  of  General  Conical  Class  have  been 
thoroughly  examined  by  investigating  mathematically  the  general  function 
as  found  by  the  Calculus  <>t  Variations. 

As  in  the  Zenithal  Class  it  has  been  found  possible  to  simplify  this 
formula  by  expansion  and  by  the  stipulation  of  making  the  total  area  of 
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the  map  true,  with  no  increment  in  the  mean  square  error  to  the  first  or 
even  second  term  of  its  expansion. 

The  Minimum  Error  Conical  Projection  of  a  /.one  is  then  found  to  be 
identical  with  the  Third  Projection  of  Murdoch  invented  in  1758. 

This  can  be  made  Orthomorphic  or  Equal  Area  as  may  be  desired, 
with  no  change  in  the  depth  of  map  or  radius  of  mid-parallel,  by  small  and 
easily  calculated  changes  in  the  radii  of  the  intermediate  parallels,  thus 
producing  two  standard  parallels  at  about  one-fifth  the  depth  of  the  map 
from  the  limiting  parallels. 

This  class  includes  the  Cylindrical  Projections,  Direct  and  Transverse, 
in  which  the  property  of  Minimum  Error  will  also  be  produced  by  making 
the  depth  of  the  map  true. 

When  the  Conical  Projections  extend  up  to  the  pole  the  principle  of 
making  the  total  area  true  also  applies,  but  the  Equidistant  Projection  is 
not  quite  so  near  the  absolute  Minimum  Error  as  in  the  zone  map.  It  is 
however  still  recommended  as  the  best  in  practice  on  account  of  its 
simplicity. 

The  general  practice  hitherto  has  been  to  make  the  maximum  error  in 
this  class  as  small  as  possible.  This  practice  leads  to  a  maximum  error 
one-twelfth  smaller  than  in  the  system  now  investigated,  but  the  average 
error  is  15  per  cent,  greater,  and  the  valuable  property  of  truth  of  radial 
scale  and  total  area  is  lost. 


The  best  projection  for  a  Cadastral  Survey. 

A  good  many  ot  the  refinements  we  have  been  considering  may  appear 
unimportant  and  unnecessary  to  the  practical  map-maker,  so  long  as  his 
work  must  be  plotted  on  a  substance  like  paper  which  changes  owing  to 
hygrometric  conditions  in  a  way  that  cannot  be  allowed  for,  from  day  to 
day,  and  even  from  hour  to  hour.  It  is  prudent  however  to  carry  our  refine- 
ments to  a  degree  further  than  the  limits  of  such  changes,  and  we  should 
always  look  forward  to  the  possibility  of  some  substance  less  changeable 
than  paper,  as  we  know  it,  being  invented. 

But  there  is  one  case  in  which  the  greatest  care  in  selecting  the  best 
projection  should  be  exercised  and  that  is  the  projection  on  which  to  base 
the  rectangular  plane  coordinates  of  a  cadastral  survey. 

In  the  great  self-governing  Colonies  and  in  most  ot  the  Crown  Colonies 
the  system  of  land  tenure  is  based  on  accurate  survey.  Boundaries  of  pro- 
perty are  defined  by  the  bearing  and  distance  of  lines  between  fiducial 
points  on  permanent  marks  set  at  intervals  round  the  property.  The  survey 
round  this  boundary  is  closed  to  a  prescribed  limit  of  error  by  the  method 
of  balancing  latitude  and  departure,  and  one  or  two  points  of  it  are  con- 
nected to  the  general  triangulation  of  the  country  either  directly  or  by 
intermediate  traverse  survey  to  a  trigonometrically  fixed  point. 

4—- 
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The  Cadastral  Survey  Department  does  not  want  to  be  worried  with 
questions  of  spherical  excess  or  excentricity  of  the  Earth,  and  in  supplying 
it  with  a  system  of  plane  coordinates  we  ought  to  ensure  that  the  errors 
of  the  projection  chosen  will  lie  well  within  the  limits  of  accidental  error 
allowed  in  the  Cadastral  Branch. 

An  illustration  from  the  Survey  on  which  the  writer  was  engaged  may 
not  be  out  of  place.  This  was  in  a  tropical  country  for  the  most  part 
undeveloped  and  covered  with  virgin  forest,  even  up  to  heights  of  5000  feet, 
consisting  of  huge  trees  from  100  to  nearly  200  feet  high,  intertwined  near 
the  ground  by  a  dense  mass  of  undergrowth  and  creepers,  through  which 
it  was  impossible  to  force  a  way  unless  preceded  by  a  man  with  a  clearing 
knife.  The  country  was  covered  with  a  network  of  major  triangles  by  the 
Trigonometrical  Branch,  of  sides  averaging  say  20  miles  in  length,  and 
very  much  of  the  interlaced  type  used  in  Great  Britain  and  Ireland.  The 
clearing  of  hill-tops  and  suitable  points  on  the  lower  ground  for  this  purpose 
involved  great  expenditure  of  money  and  labour.  This  major  triangulation 
was  however  of  very  little  direct  use  to  the  Cadastral  (or,  as  it  was  called, 
Revenue)  Branch,  as  settlement  in  such  countries  always  takes  place  at 
the  beginning  along  rivers,  roads  and  railways,  which  are  necessarily  on 
low  ground.  A  general  system  of  minor  triangulation  was  in  this  case  ruled 
out  by  the  great  cost  of  clearing.  But,  owing  to  the  development  of  the 
system  of  accurate  traversing  with  the  long  tape  of  steel  or  invar,  first  in- 
stituted in  Australia,  the  difficulty  was  overcome.  The  method  used  was 
to  run  high  class  traverses,  locally  called  standard  traverses,  along  roads, 
railways,  and  rivers  (and  sometimes  watersheds),  leaving  permanent  marks 
in  pairs  or  threes  every  half-mile  or  so.  The  traverses  were  connected  with 
the  major  triangulation  at  intervals  where  suitable  opportunities  occurred, 
generally  from  7  to  10  miles  apart,  either  directly  or  by  the  fixation  of 
suitable  minor  points.  The  limit  of  error  allowed  in  such  standard  traverses 
when  closed  on  trigonometrically  fixed  points  was  i  in  10,000,  but  the 
average  yearly  limit  attained  was  near  i  in  30,000.  These  traverses  then 
afforded  a  frame-work  ready  at  hand  in  the  most  convenient  places  for  the 
Cadastral  Branch,  with  which  to  connect  their  boundary  surveys.  The  limit 
of  error  allowed  in  the  Cadastral  Branch  was  i  in  4000,  but  the  average 
attained  was  more  like  i  in  10,000. 

The  projection  used  in  this  survey  was  the  same  as  that  used  in  Great 
Britain,  India,  South  Africa,  etc.,  namely  the  Cassini  or  Transverse  Cylin- 
drical Equidistant,  cast  on  an  average  value  of  the  Earth's  radius  in  the 
vicinity,  and  with  only  one  standard  (so-called)  parallel,  i.e.  the  meridian. 
We  have  seen  that  in  this  projection  the  .v  or  N- S  coordinates  are  always 
too  large  in  the  ratio  sec  0  :  i,  where  0  is  the  spherical  angle  of  a  point 
away  from  the  central  meridian.  When  0  is  small  the  error  is  therefore  .',6'-'. 
If  0  corresponds  with  a  distance  of  50  miles,  as  it  did  in  one  instance  in 

1  /     5°     V"  1-1 

this  Survey,  tins  error  amounts  to      I       "          I       i  in  12,300  nearly,  which 
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is  rather  uncomfortably  near  the  limits  of  error  in  the  traverse  surveys.  All 
areas  on  which  land-tax  is  paid  are  increased  to  the  same  degree,  which  is 
an  inconvenience  to  the  lessee  at  any  rate.  Of  course  this  error  was  allowed 
for  in  stating  the  closing  errors  of  the  traverse  surveys,  but  the  smaller  it 
can  be  made  the  better,  and  our  investigations  have  shewn  that  by  reducing 
all  N-S  coordinates  in  the  proportion  \  8-,  where  3  is  the  limiting  spherical 
angle  of  the  projection  from  the  central  meridian,  the  maximum  error  is 
reduced  in  the  ratio  2  :  3,  and  the  total  area  of  the  projection  at  any  rate 
is  made  correct. 

It  is  curious  that  this  device  of  making  standard  two  parallels  (or  sub- 
meridians)  in  the  Conical  or  Cylindrical  Projections  when  used  in  national 
surveys,  does  not  seem  to  have  been  noticed  before.  The  only  mention  the 
writer  can  find  of  the  use  of  such  a  device  is  by  Jordan  in  referring  to  the 
Survey  of  Mecklenburg  (p.  350)  where  the  Conical  Orthomorphic  Projection 
with  initially  one  standard  parallel  was  used,  but  by  a  little  device  (durch 
einen  kleinen  Kunstgriff),  which  was  evidently  the  use  of  a  scale  constant 

S-j 

i  —  -    ,  the  maximum  scale  error  was  reduced  one-half. 
4 

The  Cadastral  Branch,  as  stated,  was  supplied  with  the  rectangular 
coordinates  of  the  trigonometrically  fixed  points  and  the  points  fixed  in 
the  standard  traverses,  in  this  Cassini  System,  but  in  supplying  bearings 
from  triangulation  stations  such  bearings  were  given  from  a  meridian  at 

o  o  o 

right  angles  to  the  great  circle  from  such  station  to  the  principal  meridian, 
according  to  the  practice  usually  followed  hitherto.  Thus  the  bearing  from 
a  point  P  to  a  point  Q  would  differ  from  the  bearing  from  Q  to  P  by  the 
spherical  excess  of  the  space  pPQ</,  where  />  and  q  are  the  feet  of  the 
perpendiculars  from  P  and  Q  to  the  principal  meridian.  This  difference 
might  amount  to  13  seconds  of  arc  for  a  distance  PQ  20  miles  running 
parallel  to  the  principal  meridian  and  50  miles  east  or  west  of  it.  This  is  a 
small  amount  compared  with  the  accumulated  error  of  a  traverse  20  miles 
long,  but  still  such  differences  are  rather  puzzling  and  inconvenient  to  the 
Cadastral  Surveyor,  and  having  once  chosen  a  projection  it  would  be  better, 
in  the  writer's  opinion,  to  supply  all  bearings  as  well  as  coordinates,  on  the 
plane  system. 

For  a  country  nearly  round  or  square,  from  which  the  triangulation  is 
not  likely  to  be  continued  in  the  X-S  or  E-\V  direction  for  any  purpose  to 
neighbouring  countries,  the  Minimum  Error  Equidistant  Zenithal  Projec- 
tion has  many  advantages.  As  already  stated  the  maximum  error  would  be 
less  than  in  the  Cassini  in  the  ratio  8  :  3  and  the  amount  of  spherical  excess 
to  be  forcibly  dispersed  would  be  about  one-half.  The  projection  is  easily 
calculated  for  it  is  only  necessary  to  calculate  by  spherical  trigonometry 
the  bearings  and  distances  of  all  points  from  the  centre  of  projection,  and 
then  to  reduce  such  distances  in  the  ratio  i  -  }  J?  :  i.  which  means  sub- 
tracting a  constant  from  the  log  distances,  before  working  the  coordinates 
by  simple  sine  and  cosine  of  the  bearing. 
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Projection  for  Surveys  of  Hanover  and  Egypt. 

As  already  remarked  in  the  surveys  of  Hanover  and  of  Egypt  the 
Cassini  Projection  was  rejected  in  favour  of  the  Conformal,  or  Orthomorphic, 
Cylindrical  System  of  Gauss.  This  differs  from  the  Cassini  in  the  fact  that, 
in  order  to  produce  the  same  scale  error  in  the  E-W  direction  as  in  the  N-S 
direction,  the  coordinate^  in  the  former  direction  instead  of  being  the  true 
distance  on  the  sphere  is  increased  in  the  ratio 


to  0,  that  is_y  is  increased  to  y    i  +  /.-  -     ...    ,  where  A'  is  the  Earth's  radius 

in  the  vicinity.  The  projection  as  used  in  the  surveys  mentioned  has  only 
one  standard  parallel,  i.e.  the  principal  meridian.  It  could  be  improved  in 
the  same  way  as  the  Cassini  by  multiplying  both  coordinates  by  the  scale 
constant  (i  —  t!8-),  where  8  is  the  limiting  spherical  angle  from  the  principal 
meridian,  which  would  also  make  the  total  area  correct. 

In  the  Orthomorphic  Projections  the  errors  of  bearing  for  small  distances 
round  a  point  are  practically  zero  and  for  large  distances  are  much  smaller 
and  easier  to  compute  than  in  the  Equidistant  Projections.  The  errors  of 
distances  are  also  easier  to  compute.  For  these  reasons  the  writer  is 
inclined  to  think  that  for  Cadastral  Survey  purposes  these  advantages  more 
than  counterbalance  the  extra  initial  labour  of  computing  the  radial  co- 
ordinates, and  that  Surveys  about  to  start  a  new  or  revised  system  of 
projection  would  be  wrell  advised  to  consider  the  merits  of  the  Orthomorphic 
projections,  whether  direct  or  transverse,  and  to  adopt  the  scale  constant 
which  renders  the  total  area  true  and  the  average  error  a  minimum.  The 
advantages  of  the  Orthomorphic  system  have  also  become  apparent  during 
the  war  for  artillery  and  aerial  purposes,  the  new  projection  for  the  war 
maps  of  France  being  based  on  the  Conical  Orthomorphic  System  with 
two  standard  parallels,  though  with  Tissot's  criterion  of  least  maximum 
error  instead  of  the  criterion  of  least  average  error  as  recommended  in  this 
Paper.  This  map  being  more  liable  to  extension  in  the  E-AY  direction  is  on 
the  Direct  System  and  the  standard  parallels  are  actual  parallels  of  latitude. 


The  Finite  Errors  in  Conical  and  Cylindrical  Projections. 

So  far  we  have  only  discussed  the  differential  scale  errors  at  a  point, 
but  as  Minks  has  pointed  out  (Chapter  IX)  these  are  perhaps  not  of  much 
use  in  practice,  when  the  question  mostly  arises  of  finding  the  actual  errors 
in  distance  and  bearing  between  two  points  at  finite  distances  apart  on  the 
man,  as  compared  with  the  corresponding  quantities  on  the  sphere  or 
spheroid.  Formulae  for  such  finite  displacements  have  been  investigated 
to  several  terms  of  their  expansions  in  the  works  of  Jordan  and  Helmert 
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for  the  Cassini  Projection,  and  in  Jordan  for  the  Gauss  Conformal  (or 
Transverse  Cylindrical  Orthomorphic).  and  Direct  Conical  Orthomorphic 
Projections.  Jordan  also  investigates  the  coordinates  and  finite  errors  of 
the  Direct  Cylindrical  Projection,  Equidistant  and  Orthomorphic,  tangent 
along  the  prime  vertical  great  circle  at  a  point,  and  called  by  him  the 
"Querachsige  Projektion."  This,  as  he  shews,  only  differs  from  the  Conical 
Projection,  tangent  along  a  parallel  of  latitude,  by  small  quantities  depend- 
ing on  the  fourth  powers  of  the  coordinates.  The  latter  has  one  advantage 
over  the  former  in  that  the  projected  meridians  are  absolutely  straight  lines. 
In  all  these  investigations  only  one  standard  parallel  is  assumed,  but  the 
results  will  apply  to  the  case  of  two  standard  parallels  without  change  for 
the  Orthomorphic  Projections  so  far  as  the  finite  angular  corrections  are 
concerned,  for  the  scale  constant  is  applied  to  both  coordinates,  and  so  far 
as  distances  are  concerned  it  will  simply  be  necessary  to  deduct  the  scale 
constant  from  the  distance  ratio  as  given  by  the  formula  of  Jordan.  In  the 
case  of  the  Equidistant  Projections,  Direct  and  Transverse,  the  formula  for 
two  standard  parallels  will  become  more  complex  as  will  be  shewn  later, 
but  the  scale  constant  will  tend  to  reduce  the  amounts  of  the  corrections. 

A  mathematical  paradox  in  the  General  Formula  of 
Minimum  Error  Conical  Projections. 

There  is  a  rather  curious  mathematical  paradox  in  the  general  formula 
for  the  Minimum  Error  Conical  Projections  which  deserves  notice  for  its 
mathematical  interest,  though  it  is  not  a  question  that  need  trouble  the 
practical  map-maker.  It  will  have  been  noticed  that  in  the  solution  of  the 
general  case  of  a  map  covering  a  /one  from  a  lower  limit  <t  to  an  upper 
limit  (3,  the  Calculus  of  Variations  exacts  that  the  radial  scale  error  shall 
vanish  at  both  limits,  i.e.  dr\dQ  =  i  both  when  0  a  and  61  ft.  If  however 
we  make  a  -  o  we  approach  the  case  of  map  extending  to  centre  and  in  the 
case  of  the  Zenithal  Projection,  as  amended  by  Clarke,  we  get  for  the  radial 
scale  at  centre,  not  the  value  i,  but  the  same  as  the  perpendicular  scale 


Again,  in  the  General  Conical  Class  the  condition  of  Minimum  Error 

rV 

leads  to  two  standard  parallels  and  a  mean  square  error  =  :   but  it  we 

\  45 

start  with  the  assumption  that  the  map  extends  to  the  centre,  we  are  led  to 
a  projection   with   one    standard    parallel   only   and   a    mean    square    error 

-f—  -       Now   if  in   the   former  value  for  /-'.  we  allow  a  to  approach  the 


value  zero,  then  f>  approaches  the  value  $ft  and         _  approaches  the  value 

s  4o 

J!L  _    which  is  slightly  larger.     Professor  Nicholson   was   kind  enough  to 
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give  some  attention  to  this  paradox,  and  without  going  into  a  full  investi- 
gation, the.  explanation  seems  to  be  that  so  long  as  we  have  a  value  for  a, 
however  small,  the  projection  will  have  the  radial  scale  true  at  a  and  two 
standard  parallels  ;  but  when  a  is  actually  zero  a  discontinuity  takes  place 
and  we  have  a  projection  with  only  one  standard  parallel.  The  paradox 
seems  to  be  somewhat  similar  to  that  of  the  deflection  of  a  circular  elastic 
plate  supported  round  the  edge.  If  there  is  a  hole  in  the  centre  however 
small  we  get  a  different  result  from  that  derived  from  the  initial  assumption 
of  plate  continuous  to  the  centre. 


CHAPTER  III 

THE  CONICAL  ORTHOMOKPHIC  PROJECTION  WITH  TWO 
STANDARD  PARALLELS  CLAM  HERTS  SECOND)  EOR  THE 
SPHEROID 

We  have  hitherto  considered  the  Earth  to  he  a  sphere  of  radius  A\  taken 
as  unity  in  our  formulae  for  brevity,  but  us  it  is  usual,  especially  in  the 
projections  employed  for  the  rectangular  coordinates  of  triangulation,  to 
take  account  of  the  excentricity,  \ve  will  now  investigate  the  expansions  of 
the  spheroidal  formulae  for  the  above.  The  expansions  of  the  formulae  for 
the  case  of  one  standard  parallel  are  given  by  Jordan  (pp.  521-4),  Pixxetti 
(pp.  389-91)  and  Oscar  S.  Adams,  Special  Publication  No.  47  U.  S.  Coast 
and  Geodetic  Survey. 

The  formula  for  the  spheroid  in  terms  of  our  original  notation  is 

tie 
,  /I  +  CCOS  0\  -2 

r  =  2 A  (tan  !  &)«  , 

\  i  —e  cos  v/ 

where  e  is  the  Earth's  excentricity.    It  will  be  now  more  convenient  to  put 
.AT  for  the  constant  2 A  and  write 


Calling  the  radius  of  curvature  to  the  meridian  p,  and  the  grand  normal  or 
radius  of  curvature  at  right  angles  to  the  meridian  r,  we  have 

p  =  a  ( i  —  e'-}:  ( i  -  e-  cos"  0) :-; ,    and    v  =  a  (i  —  r  cos'-  6>y:', 

where  a  is  the  semi-axis  major  of  the  spheroid.  The  scale  along  the  meridian 
or  what  we  have  hitherto  called  the  radial  scale,  is  given  by  dr  dm,  where 
dm  is  an  element  of  the  projected  meridian.  Now 

i  dr      dr  (i  —  r  cos'J  &)- 

dr  dm  =  drdO.dQ  dm  =  -  .,'. 

/j  dti      do       a  (i  —r) 

Differentiating  logarithmically  we  have,  putting  r'  for  dr  dt),  >-"for  d~r  d&~,  etc., 
log  /-  =  log  J\  +  a  log  tan  \  6  +  J  [log  ( i  +  e  cos  0)  -  log  (i-e  cos  0)\, 


Therefore 

and  as  the  projection  is  orthomorphic  the  scale  along  the  parallels  has  the 
same  value.    Putting  as  before  x  for  the  co-latitude  of  the  mid-parallel  and 
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8  for  the  angle  from  x  to  standard  parallels  0n  02,  we  have  0i=x~^ 
02  =  X  +  g.  The  condition  that  the  scale  at  Ol  shall  be  the  same  as  the  scale 
at  00.  is 

( i  -  ?  cos-  0,  )4  ( i  -  r  cos-  02)i 


ri.f  -i  >          /i  —   "•'  •>  /i 

(7  sin  0j  a  sin  02 

where  r1;  ;-.,  are  the  radii  at  0]}  0,.    Writing  out  the  full  expressions  for 
rlt  r.2  and  cancelling  common  factors  this  leads  to 

log  sin  0o  -  log  sin  Ol  -  log  ( i  -  e~  cos-  0.,)^  +  log  ( i  -  e-  cos'-  0,)^ 

+  t'cos  0A-2 


The  additional  condition  that  the  scale  shall  have  the  value  unity  along 
these  parallels  is 


which  leads  to 

v,  sin  0,      /i  —  t'cos  0AT  p.,  sin  0.,      / 1  —  c  cos  0.A  T 

A  —  — 


(tan  !  61!)"  \i  +  c  cos  #,/          //  (tan  '-  O,)'4  \i  +  6'  cos 

It  will  also  be  noted  that  the  radii  of  the  standard  parallels  have  the  simple 
values  r1--=v1sm61/Jt,  r.>  =  r.2  sin  O.,jn.  Such  are  the  complete  formulae  for 
two  standard  parallels  chosen  beforehand.  To  calculate  the  values  of  r  by 
them  with  exactitude  log  tables  of  at  least  S  figures  are  required,  so  that  in 
practice  it  is  more  convenient  and  even  more  accurate  to  use  expansions, 
involving  the  differences  and  values  at  the  mid-parallel.  These  expansions 
can  be  made  in  terms  of  either  the  difference  of  latitude  A0  following 
Jordan  or  the  difference  of  meridian  distance  Aw  following  Pix/etti,  the 
U.S.  Survey  and  others.  These  writers  consider  the  case  of  one  standard 
parallel  only,  where  n  =  cos\, 


//  (tan  .',  x)"  V  i  -r  e  cos  x 
It  may  be  remarked  that  if  we  write 


then  £  is  very  nearly  the  geocentric  co-latitude 

=  6+  .Irsin  20-:-  J-^sin  26  cos-  0+  ..., 
whereas  by  expanding  the  above  expression  we  find 

£  =  6  -i-  !  e"  sin  28  +  ,"'.,  c1  sin  26  cos-  6+  ... 

as  pointed  out  by  Her/  (p.  246).  If  tables  of  geocentric  co-latitude  are 
used  therefore  which  take  account  of  the  term  in  c4  the  error  is  very  small,, 
the  maximum  being  only  o"'5  for  0  30°. 
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By  expanding  the  expressions  involving  e  we  can  write 
_  Iog10sin#2-log,0sin  0J 


~  loglt)  tan  £  02  -Togu,  tan  -J,  0^1/er  (cos  0.]  -  cos  0,)  +  .'^/^(cos^-cos^,  )+...' 

which  is  fairly  easy  to  compute,   M  being  the  modulus  of  common   log- 
arithms. 

But  for  our  purpose  we  require  to  expand  the  whole  expression  for  n  in 
terms  of  8'-'  and  e~  which  can  be  done  as  follows  : 


,-A, 

log  sin  0o  —  log  sin  6^  ~  I     cot  6dO 

Je, 


«  /  8-        I  84     2  +  COS- 

:2Scotx  (  i  +-—-—+ 

V       3  sin-x      15      sin4x 

/•ftj 
log  tan  -',-  0.,  -  log  tan  !  6}  =     '  cosec  Bdtt 


/       8-  T  -<-  cos-  v       84    ?  +  1  8  cos-  Y  +  cos4  v 
=  28  cosec  Y    i  +  f  ^  +  - 

\        0      sm'x          i  20  sin  ^ 


i    -  I  2  0 

ir(cos'-0o-cos-01)=-r  I     sin0cos0d?0  =—  2oe~  (  sinxcosx-   --  sinx  cosx 
h.  3 


2S2    . 


We  shall  not  carry  the  expansions  beyond  84  and  e~82,  but  if  desired 
the  further  proceedings  are  obvious.     Dividing  we  find 

2  0  COt 


I        6  v  360         sin-x 

WTe  shall  now  continue  the  expansion  of  /-,  but  it  will  be  convenient  to 
do  this  in  terms  of  Aw,  the  meridian  distance,  instead  of  A0,  the  difference 
of  co-latitude.  We  have  found  already  dr[dm  ur(\  e1  cos"  0)-  a  sin  0, 
and  differentiating  further  logarithmically,  dropping  the  factor  a  pro  te/n., 

we  have 

r"      f  r  c  -  n  cos  0  ( i  -  e'-) 

(-   I  '  -i  ' 

^       \r/  sin-0(i      t'-cos-^)- 

and  multiplying  by  d&  dm, 

/r"      (r'\'\         ~  "  cos  ^  ^z     e~  COS'J  ^' 
//-i  sin'J^ 

/  f"\  !)    Ill    —    C(  IS    D  I    I 

therefore 
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Continuing 

r"      r"  r          (i  —  e"  cos2  0      (i  — <r)  2cos0(«  — cos0) 


r         r1  I       sin  0  sin3  0 

r"      r" r'\  f(i  —  tr cos20):^      2  cos  0  (n  —  cos  0)  (i  —  r 


But 
therefore 


r         r8  /»,   '     \sm0(i-<?)  sin:i0 

y>'\         n-  (n  -  cos  0)  ( i  —  f  cos2  0)- 


Again 


r    r          ($e-  (i—e-  cos'-'  6)*  cos  0  sin'2  9  —  cos  0  (i  -  e"  cos2  0)- 


/•          r'2  1  sin"  0  (i  —  t'- 


y'2  cos  0  (i  —  <r  cos2  0)-  (//  —  cos  0)  (n  —  2  cos  0) 

sin2  0 
3  cos  0  (//  —  cos  0)  (n  —  2  cos  0)  ( i  —  £2  cos2  0)^ 

and 
riv      /'V\  j-  cos  0  (i  —  r  cos2  0)4  +  5r  cos  0  sin2  0  (i  —  f-  cos2  0):! 


sin2  0(i  -t-2)2 

(/7  —  2  COS  0)  (i  -  C2  COS2  0):!         2  (/V  -  COS  0)  ( I  —  C1  COS2 

sin2  0  ( i  -  e'2)  sin2  0  ( i  —  ^2) 

3^2  cos  0  ( i  —  e~  cos2  0)2  (n  —  cos  0)  (11  —  2  cos  0) 
sin2  0  (i  —  e1) 

3  cos  0  (n  —  cos  0)  (;/  —  2  cos  0)  (i  —  r'  cos2  0):i 
Now 


sin'-  u  (i  —<•'-)- 

/  Zj\    /  "  •'  Zj \ " 

4  (n  —  cos  (7)  ( i  --  e~  cos^  r/)" 

sin2  #  ( i  -  r'J) 
3t'2  cos  0  (;/     cos  0)  (//    -  2  cos  0)  ( r  —  r1  cos'-'  0):; 

sin2  0  (i  —  6'-) 

3  cos  0  (;/  -  cos  0)  (//  -  2  cos  0)  ( i  —  e1  cos2  0):; 
~  sin4  0(i-  7)~ 
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—  cos  0  (  i  -  e-  cos2  0)'-'  (  i  —  e'2  -  4^-  sin-  0} 

sina  0  (  i  -  <?-)- 
4  (n  —  cos  6*)  (  i  -  e-  cos'J  0):i 

si  11*0(1  -<T~ 

3  cos  6  (>i  -  cos  0)  (/i  -  2  cos  6)  (  i  -  <~  cos-  0)- 
sin4  6 

n(n-  cos  0)  («  -  2  cos  0)  (  i  -  f"  cos'-'  6)-  ) 
sin4  ft  /  ' 

Now  if  we  are  expanding  only  as  far  as  m4  we  can  put  n  cos  \  in  the 
expressions  for  (r"  jr}m  and  (r'^r),,,,  so  that  only  the  first  terms  within  the 
brackets  will  be  retained.  Since  also 


we  have 

pv  =  a'-  (  i  —  6'-)/(  i  -  t'-  cos-  fr)'J,     pi/-  --  -  a'''  (  i  —  fc1-),  (  i  -  r  cos-  0)  -. 

Thus,  replacing  0  by  the  value  x  for  mid-parallel,  we  have 


( i  -  e~  cos-  \-    i  --  e-  -  ^r  siv  x     >«       ) 
a4sin'Jx  (i  -fL')-  24      J 

We  shall  subsequently  expand  K  in  terms  of  x  and  8  but  for  the  moment 

we  put 

ne 

v,  sin  0.,      /r  —  e  cos  #./,  - 
" 


We  then  have  radius  at  mid-parallel 

a  sin  x 
r  -=  ,  (i  -:-  terms  in  o), 

H(I  -rcos2x)a 

and 

|-         (r  -  cos  x)  ( i  -  ^cos2  x)--    nr      ( \  -  r  cos'-  \f   m:' 


_co^(i_^rcos^/    _  4£^xV"       1        +  icfms  m  g) 
(?:;(l  -6'-;  i  -^    /  ?4 

|"         (;/  —  cos  x)     M~  _^   l     m" 
^  (_  i/sinx      '    2       pv    6 

cot  Y  /        4<J  s'nJ  \\  »'4       \    •  • 

_    A  (  i  -  4  A          ...  V  ( r  +  terms  in  o). 

pv2    V  i  - ''"    /  24      i   . 
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For  n  -  cos  x  this  reduces  to  the  expression  found  by  Pizzetti  for  one 
standard  parallel,  when  terms  in  8  vanish.  Substituting  however  the 
approximate  value  of 


we  find  the  second  term  gives  a  term  in  w2S-  and 
r-rx  =  {«+-.£  +  ^  (i  +  ^sin'x)  ^ 

COt  X  ,  2    •    2     N  "l*         \  1 

-,-  ( i  —  ±e-  sin-  x)  —  •  •  •  M  i  +  terms  in  o). 
Pv  24       ) 

We  must  now  expand  the  expression  for  K  in  order  to  find  the  terms 
in  S.    We  have 

ne 

a  sin  6>  (cot  i  ft)"  / 1  -  £  cos  ft,\  a 
K^ 


Now  sin  0.,  =  sin  x  (  J  +  cot  X  •  ^  ~  ~  +  T  cot  X  •  •  • )  > 

also 

/          x/     „     ;/-  +  n  cos  x  <$- 
(cot  -?,  ft,)"  =  (cot  ?,  x)    (  i  -  -r—  8  +  — r-5 .  - 

V       sin  x  sin-  x         2 


and  putting  in  the  proper  value  of  n  this  becomes 
(cot  \  ft,)"  =  (cot  I  x)"  (  i  -  cot  x  •  S  +  cot"x  •  s" 

-C0tx(2-j^f: 

Multiplying  this  by  the  expansion  for  sin  ft,  we  find 

sin  ft,  (cot  \  ft,)"  =  sin  x  (cot  \  x)"'  f  I  -     -  'f-  r  sin  x  cos 
Next  take  (i  —  e-  cos'J  ft,)    a  =  y  say,  we  have 


order,  and      -  =  ^c'-  sin  61  cos  61  also  to  same  order,  therefore 
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In  the  same  way  we  shall  find 


i  -  e  cos  (9.A  -        /  1  -  e  cos  Y\  - 

- 


. 
i  -r  e  cos  EV          \i  +  e  cos  \ 


\  -   /  tf 

)     I  I  4-  nc-  sin  vd  +  nr  cos  Y  .  -- 
l     \  *     2 


/  1  -  ^  cos  x\  -  /  8-  \ 

"  \  i  +  fcosx/     \  l  +  e'  SU1  x  COS  x  •     +  *'"  COS"  x  •  -,  +  terms  in  ^  upwards  J  . 

Multiplying  out  we  get 

nc 
/I-gCOsgB\3 

—  ^      (i  -  £-  cos2  ft,)  ~  2 

\I  4-cCOS  6J 

»(' 

(\  -ecosx\'2  ,         .,       .,    .     ,  /  3"  28"         \ 

(  i  -  r  cos-  Y)     2(1  +  r  sur  Y  .  -  +  r  sin  Y  cos  Y  .        4-  ...  I  , 

\I   -rd'COS  X/  V  2  A       3  /  ' 

and  multiplying  this  by  the  expression  above  we  get 


^  -  ism 

«(i-rcos-X)"  ^  +*COSX/ 

+  terms  in  84  upwards]. 

We  can  test  the  accuracy  of  these  expansions  by  seeing  if  the  scale 
works  out  correct  on  the  standard  parallels  ;  we  have 


dm  v  (       Pv 

cot 


(i  —  r  cos2  xV-  (i  -  4?"  sin'-  Y)  cot  x   »i''       }   I        S-  ) 
^—r-, -, -•  —  ...  i    LI-  -0  -6-siirx)...f. 

Now 


Substituting  these  and  neglecting  powers  of  c  above  c-  we  get  for  the  scale 

T  x 5"  ) 

J  i  +  ( r  -  t"  sin-  x)      -  <~  sin  X  co-s  X  •  6'  *  •• 

I  2 

which  proves  our  formulae. 
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If  we  write  m0  for  the  mean  of  the  meridional  distances  from  mid-parallel 
to  standard  parallels  we   have   m0  =  p8  +  terms   in  r'S:i   upwards,  so   that 

S2  =  —  §  to  same  order  and 
P 

8»  (i  -  *  sin-'  x)  =  <ilrr2  sin  *)  (i  -  ^  cos'  X)'  =  < 

a-  (  i  -  t'-)-  pv 

up  to  terms  in  e*  inclusive,  so  that  we  can  write 

»?         n 
+  ?  - 
6pv      24p 


{,-*•'...}. 

(.  2pv          J 


The  value  of  rK  is 


< 
=  v  tan  x  j  i  -  —  (i  +  \r  sm- 

V.  O 

The  French  in  using  these  expansions  for  their  new  maps  employ  up 
to  the  term  mAj6pv  only,  and  for  this  purpose  it  is  sufficient  to  take  for  n 
the  value  cos  x  as  for  one  standard  parallel,  the  correction  to  make  two 
parallels  standard  being  effected  by  the  multiplier  (i  -  ;//02/2/>v).  For  r^ 
they  use  the  value  rtanx(i  —  w(,-/2pv).  By  some  oversight  they  employ  v1 
instead  of  pv  in  their  formulae.  These  small  errors  by  a  fortunate  con- 
currence of  opposing  effects  have  an  almost  imperceptible  result  on  their 
tables  of  coordinates. 

The  writer  has  expanded  the  formulae  still  further  as  for  the  sphere, 
that  is,  neglecting  the  terms  depending  on  e.  For  r—rv  there  should  be 


added  +  -  7-r,,  ----- ,  and  the  reducing  factor  becomes 

I20A4     sm'x         36A>4sm-x 

rr  .  84  i  -  2  cos"  Y 

i  -   -  ( i  -  c-  sin'-  v)  +  —'•')  tne  further  term  in  n  has  already 

2  v  24        sin'x 


been  given.  If  we  take  x  :  9°°  an^  Put  smX==I»  cosx-o  we  obtain  the 
formulae  for  the  Mercator  Projection  for  the  spheroid  with  two  standard 
parallels  in  latitudes  +  8. 

Minimum  Error  formulae  for  the  Spheroid. 

In  the  case  of  the  minimum  error  formulae  for  the  spheroid  we  do  not 
select  the  two  standard  parallels  beforehand,  but  we  make  the  scale  values 
equal  on  the  upper  and  lower  parallels  for  one  condition  determining  ;/; 
and  for  the  other  determining  A' we  make  the  total  depth  of  the  map  true. 
Now  we  have  already  in  the  above  investigation  determined  the  formula 
for  n,  because  we  have  only  to  consider  8  as  being  no  longer  the  angle 
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from  mid-parallel  to  standard  parallels  but  as  the  angle  from  middle  to 
extreme  parallels.    To  find  K  we  have 

ne  in 

+  e  cos  8\  2  .„  /i  i-  c  cos  u\  -  )        ft 

--£)    -(tanJa)H(--  =       pdB. 

—  ccos/v  \i-ecos  a/    J       Ja 

We  have  already  indicated  the  method  for  expanding  these  expressions, 
which  involve  long  and  complicated  differentiations  and  it  will  suffice  here- 

to give  the  result, 

ne 

^     v  sin  x  (cot  .|x)"  /i  -  e  cos  XV  f  i  _  m,':       j 
n  \  i  +  c  cos  x/     I         6pi/      J  ' 


so  that  our  reducing  factor  is  now  (  i  —  1     -  )    where  ;;/„  is  the  mean  of  the 

V  pvj 

meridional  distances  of  upper  and  lower  parallels  from  the  mid-parallel. 
The  corresponding  value  of  rx  is 

a  sm  X          f        8"  /  1  f        8-         .....          ] 

-T    i  -  :r(i~rsmLx)---[  =vtanX  |i  -  -  (i  +  :J  -sm-x)  •••     , 
«  (i  —  rcos2xr  3 

where  8  is  now  the  angle  from  middle  to  extreme  parallels. 
It  is  obvious  that  by  writing 

n?  (        IH{?  } 

r  -  ;-  =  ;;/  —          and    A  =  ...  li  +  -f  —  ,-  , 
opv  (         opi'J 

we  shall  obtain  a  projection  which  is  (approximately)  Equal  Area  such  as 
would  result  from  expanding  Albers'  minimum  error  formulae  for  the 
spheroid. 


CHAPTER  IV 

THE  POLYCONIC  PROJECTIONS 

The  Simple  Poly  conic  Projection,  which  the  French  cartographers  call 
the  "projection  polyconique  ordinaire  des  Americains"  because  it  seems 
to  have  been  first  introduced  by  the  U.S.  Coast  and  Geodetic  Survey,  has 
the  following  characteristics.  A  standard  meridian  of  longitude  is  chosen 
which  is  usually  the  centre  line  of  the  map.  This  is  developed  in  its  true 
length  on  the  map  and  each  parallel  of  latitude  is  developed  on  its  own 
tangent  cone  opened  out  on  the  plane  of  the  map.  Thus  the  standard 
meridian  is  true  to  scale  on  the  map,  and  all  the  parallels  cut  it  at  right 
angles  and  preserve  their  true  scale ;  but  because  the  radii  of  the  projected 
parallels  vary  as  the  cotangent  of  the  latitude  from  point  to  point,  the 
developed  meridians  on  each  side  of  the  standard  are  no  longer  straight 
but  curved  lines  concave  to  the  central  meridian  and  their  lengths  are  too 
great  in  proportion  as  their  difference  of  longitude  from  the  central  meridian 
increases.  When  this  projection  is  used,  however,  for  only  a  small  distance 
on  each  side  of  the  central  meridian  its  errors  are  correspondingly  small, 
and  it  has  the  advantage  that  concise  Tables  can  be  calculated  once  for  all 
for  plotting  the  coordinates  of  the  graticule  of  any  map.  Extensive  Tables 
of  this  kind  have  been  published  by  the  U.S.  Coast  and  Geodetic  Survey, 
for  each  minute  of  latitude  from  the  Equator  to  the  Pole,  and  for  intervals 
of  longitude  up  to  30°  from  the  central  meridian,  in  the  Report  of  that 
Survey  for  1884.  Full  particulars  of  the  theory  as  applied  to  the  spheroid 
will  be  found  in  Special  Publication  No.  57  of  that  Survey.  These  Tables 
were  not  used  by  the  Committee  for  the  International  Map,  as  different 
values  of  the  Earth's  spheroidal  elements  were  thought  desirable.  It  seems 
to  the  writer  very  desirable  that  some  international  agreement  should  be 
come  to  in  matters  of  this  kind,  as  many  valuable  Tables  have  been 
calculated  in  different  countries  using  different  values  of  these  elements, 
thus  seriously  vitiating  their  usefulness  for  purposes  of  comparison.  As  the 
actual  geoid  differs  considerably  in  places  from  any  of  the  spheroids  put 
forward  from  time  to  time,  but  as  these  spheroids  differ  among  themselves 
by  quantities  so  small  as  hardly  to  concern  the  map-maker,  the  best  course 
would  appear  to  be  to  decide  on  which  of  them  the  most  work  has  already 
been  clone  in  the  way  of  calculating  Tables  and  to  adopt  that  spheroid  for 
all  international  purposes. 

The  equivalence  of  the  Simple  Polyconic  Projection  with  the  Cassini 
will  be  readily  seen  by  comparing  the  expressions  for  the  expansions  of  the 
coordinates.  Those  for  the  Cassini  art:  given  by  Jordan  (p.  483).  The 
origin  is  supposed  to  be  at  the  Equator,  .v  is  the  coordinate  in  the  NS  and 
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y  the  coordinate  in  the  E-W'  direction  positive  N  and  E,  9  is  the  latitude 
and  X  the  longitude  (reckoned  from  the  central  meridian  in  c.m.)  of  a  point 
on  the  map,  M  is  the  meridian  distance  of  the  point  from  the  Equator 
(called  B  by  Jordan),  v  the  grand  normal  at  latitude  9.  Then  we  have 

Cassini       x  =  Af+  \,  v\-  sin  9  cos  9  +  ./t  i'X4  sin  9  cos  9(5-6  sin-  9)1-  . . . , 
Polyconic  .v  =  J/+  \  rX-  sin  9  cos  9  -  ^  i/X4  sin:!  9  cos  9  . . ., 

Difference  ./',-vX4  sin  <£  (-Os::  4> 

Cassini  y  =  rX  cos  9  -  ,1  i/A:!sin-<£  cos  9  -  ,  ] ,-,  rXn  sin-  9  cos  9(8-9  sin'9). . ., 
Polyconic  y  =  vX  cos  9  -  ,lrX:isin-ri  cos  9  +  , .',  ,-,-i/X5  sin'  9  cos  9  ..., 

1  )ifference  ,'.  rX"'  sin-  9  cos:!  9  .... 

In  the  International  Map  of  the  World  on  the  scale  of  i. 'million  the 
device  was  adopted  of  reducing  the  scale  along  the  central  meridian  by 
half  the  amount  of  the  error  at  the  extreme  meridians,  that  is,  by 
-|  X- cos'-' 9,,, ,  where  9,,,  is  the  mean  latitude  of  a  sheet.  This  appears  to 
have  been  regarded  as  a  rather  empirical  correction  by  its  proposers,  but 
the  investigations  we  have  made  prove  that  it  was  really  based  on  a  logical 
criterion,  that  of  making  the  maximum  error  as  small  as  possible.  In  the 
writer's  opinion,  however,  it  would  have  been  better  to  take  as  the  reducing 
factor  ,'.  X- cos'-' 9//( ,  following  Airy's  criterion  of  making  the  average  error  a 
minimum  and  the  total  area  true.  But  as  a  matter  of  fact  the  correction  is 
on  these  sheets  so  small  as  to  be  within  the  errors  caused  by  hygrometric 
changes  in  the  paper.  It  may  be  noted  here  that  in  these  Instructions  and 
in  Germain,  .v  is  used  for  the  E-W  and  y  for  the  N— S  coordinates.  The 
writer  prefers  to  use  these  symbols  the  other  way  about,  for  if  a  is  the 
bearing  of  a  line  s  nearly  always  reckoned  from  the  North  in  survey  work, 
we  have  the  common  forms  .v  =  s  cos  a,  r  =  .v  sin  a. 

Equal  Area  and  (approx.)  Orthomorphic  Polyconic  Projections. 

T.  I.  Craig  has  given  on  p.  6r  of  his  Theory  of  J\fap  Projections  the 
formulae  for  a  Polyconic  Projection  which  is  Equal  Area.  Instead  of  the 
inclination  of  the  projected  meridians  being  X  sin  9  as  in  the  Simple 
Polvconic,  this  quantity  becomes  a  series  of  which  the  first  two  terms  are 

/         rX-  cos'-  9       \ 
X  sin  (fill  - 

\  op  / 

By  changing  the  sign  of  the  second  term  and  taking  this  quantity  as 

X  sin  4>  ( i  '    •••  )   Wd  obtain   a  polyconic  projection  which  is  ap- 

6p 

proximately  orthomorphic.  It  will  be  found  to  agree  in  fact  with  the 
expansion  of  the  Transverse  Mercator  or  Gauss  Contormal  Projection  to 
the  first  term  of  the  scale  errors. 

It  would  appear  to  be  impossible  to  obtain  a  truly  Orthomorphic 
Polyconic  Projection  without  sacrificing  uniformity  of  scale  along  the 
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central  meridian.  Such  a  projection  was  invented  by  Lagrange  and  is 
usually  known  by  his  name.  That  it  really  belongs  to  the  Polyconic  Class 
is  pointed  out  by  Tissot  and  Herz,  and  also  by  Oscar  S.  Adams  in  Special 
Publication  No.  57  of  the  U.S.  Coast  and  Geodetic  Survey. 

Lagrange's  Projection  for  the  Spheroid. 

The  full  investigation  of  Lagrange's  Projection  for  the  spheroid  will  be 
found  in  Fiorini  (pp.  533  et  seq.)  and  U.S.  Publication  No.  57,  and  it  is 
not  necessary  to  repeat  it  here.  The  principal  meridian  coordinate  x  can 

A(kc"  -tn\lk 
be  written  in  the  form  .r  =  ,  ,,  ',,  ,  where  A.  k  and  n  are  constants 

kcn  +  2  +  tn[k 

lie 

to  be  determined  and  cn  -  i/f"  is  written  for  (cot  \nY  (  -~nu  being 

\i 


the  co-latitude  of  a  point  on  the  meridian.  The  scale  value  dxjdm  will  be 
given  by  2>iAtv  sin  u  (kcn  +  2  +  t"jk),  and  Lagrange  chose  the  constants  k 
and  n  so  that  the  first  and  second  differential  coefficients  of  this  quantity 
should  vanish  at  the  point  of  co-latitude  ?/„  chosen  as  the  centre  of  the  map. 

...  n  +  cos  UH 

I  he  first  condition  leads  to  KC™  =  .    From  the  second  condition 

n  —  cos  ?/„ 

Lagrange  derives  n-(i  +  sin2?/,,)"-,  but  in  doing  so  he  neglects  the  excen- 
tricity.  Captain  G.  T.  McCaw  and  the  writer  almost  simultaneously  found 
the  correct  value  for  n  should  be  n  =  |i  +  sin2  ?/0  +  2e~  sin4  ?/0/(i  -  e'-)fi.  The 
constant  A  is  usually  chosen  so  that  the  scale  shall  be  true  at  co-latitude  ?/„, 
but  our  investigations  shew  that  a  better  result  would  be  achieved  by  taking 
A  so  that  the  scale  error  along  the  central  meridian  should  be  one-third 
that  of  opposite  sign  along  the  extreme  meridians. 

For  a  map  of  small  extent  in  longitude  Lagrange's  formulae,  when 
expanded,  agree  to  the  first  term  of  scale  errors  with  those  of  the  Transverse 
Mercator,  so  there  is  no  advantage  in  using  it  for  maps  of  moderate  size. 
The  variation  of  scale  along  the  central  meridian  is  very  slow.  The  writer 
finds 

A?/3      3  —  2  sin-  ?/,,     A//4       \ 

A.Y  =  A  A//     i  -•-  cot  //„  -  -  .    ,  .  •  •  •     , 

\  24  sin-  ?/„          1 20       / 

so  that  if  the  centre  of  the  map  is  on  the  Equator  even  the  third  differential 
coefficient  of  the  scale  value  vanishes. 

This  projection  has  some  interesting  geometrical  properties,  as  the 
projected  meridians  and  parallels  are  true  circles  even  on  the  spheroid.  It 
has  this  advantage  over  the  Mercator  that  a  map  of  the  globe  can  be 
constructed  on  it  without  the  scale  error  becoming  infinite  (which  applies 
also  to  all  the  Polyconic  Projections)  but  the  exaggeration  from  80°  latitude 
to  the  Pole  is  enormous  (compare  Projections  Xos.  v  and  Yin  in  Germain, 
Plates  VI  and  VII;. 


CHAPTER  V 

THE   FINITE   ERRORS   OF   PROJECTIONS 

As  we  have  already  remarked  in  ('hap.  II,  the  differential  scale  errors  at 
a  point  of  a  projection  are  not  perhaps  of  much  value  to  the  users  of  maps. 
Generally  speaking,  the  projection  should  be  so  chosen  that  the  greatest 
error  in  distance  and  bearing  between  two  points  at  a  finite  distance  apart 
taken  from  the  map  shall  be  smaller  than  the  accidental  errors  due  to 
stretching  of  the  paper,  etc.,  and  smaller  than  the  degree  of  accuracy 
required  in  the  purpose  for  which  the  measurement  is  required.  Hut  when 
the  projection  is  required  for  the  purpose  of  fixing  points  svith  great  accuracy 
by  rectangular  coordinates  as  for  a  Cadastral  Survey,  or  even  for  artillery 
work,  it  may  be  necessary  to  know  how  the  distance  and  bearing  derived 
from  such  rectangular  coordinates  differ  from  their  true:  values  on  the 
surface  of  the  sphere  or  spheroid. 

l;or  the  Gauss  Conformal  or  Transverse  Mercator  Projection,  as  used 
in  the  Surveys  of  Hanover  and  Egypt,  such  formulae  have  been  worked  out 
for  the  spheroid  to  a  high  degree  of  accuracy  by  General  Oscar  Schreiber 
in  his  TJicoric  der  Projektiousmethode  dcr  Hannoverschen  Landesvermessung 
(Hannover.  HahiVsche  Hofbuchbandlung,  1866).  They  are  also  given  as 
far  as  the  fourth  powers  of  coordinates  by  Jordan  (pp.  507-12).  These 
formulae  apply  to  the  case  of  scale  true  along  the  meridian  (or  what  we  have 
called  one  standard  parallel  in  our  investigation  of  the  Conical  Projections), 
but  in  the  case  of  orthomorphism  they  can  be  applied  to  the  similar  projection 
with  two  standard  parallels  or  sub-meridians  without  any  change  in  com- 
puting the  angular  corrections,  and  for  the  distance-  correction  it  will  only 
be  necessary,  after  computing  the  correction  by  these  formulae,  to  apply 
with  opposite  sign  the  scale  reduction  factor. 

The  same  formulae  can  be  used,  with  simple  interchange  of  the  .x-  co- 
ordinates for  the  r,  at  any  rate  to  the  first  term  of  the  errors-— which  ought 
to  be  near  enough  for  all  practical  purposes- -to  the  Direct  Mercator.  the 
Direct  Conical  or  Lambert's  Second,  anil  the  ( 'ylindrical  Tangent  to  a  Prime 
Vertical  (Jordan's  Ouerachsige)  Projections. 

As  to  the  Cassini,  standard  along  one  meridian.  Jordan  besides  others 
gives  the  formulae  for  the  finite  errors  for  the  sphere  but  not  for  the 
spheroid.  These  formulae,  to  the  first  term  at  any  rate,  can  tie  applied  to 
the  spheroid  by  substituting  i>r  for  K',  where  />  and  i  are  taken  tor  the  mean 
latitude  of  the  two  points  in  question.  The  correct  tormulae  tor  the 
spheroid,  as  far  as  the  fourth  powers  of  the  coordinates,  are  given  by 
Helmert  (equation^  (5)  and  (6),  p.  417),  after  a  rather  difficult  investigation 


;o 
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and  in  a  rather  difficult  form.    The  writer  succeeded  in  checking  Helrnert's 
result  in  two  ways,  first  by  a  transformation  from  the  coordinates  of  the 


Fig.  2. 

Gauss  Conformal,  and  secondly  by  pursuing  the  method  of  Gauss  (see 
Schreiber,  p.  39)  of  forming  the  differential  equation  of  the  curve  assumed 
by  the  projected  geodesic  between  the  t\vo  points,  by  the  Calculus  of 
Variations.  This  differential  equation  is  more  complex  than  that  in  the 
Conformal  System,  but  a  solution  can  be  found  in  series.  It  can  be 
written  thus: 


the  symbols  having  the  same  meaning  as  in  Jordan,  p.  511,  where  an 
abridgement  of  Schreiber's  method  is  given.  This  equation  when  solved 
will  give  the  angles  between  the  chord  joining  the  two  points  on  the 
projection  and  the  tangents  to  the  projected  geodesic  ;  but  to  these  must  in 
addition  be  applied  the  change  in  the  bearing  of  this  chord  due  to  the 

projection  not  being  orthomorphic  =  ,  A^vAr.     Finally  we  shall  arrive 

"*  7"  -  s'- 

at  the  following  equations,  where 

A.v  =  .v._,  -  .v, ,      Av  -j>2  —  j'j ,      s~  —  AA-  +  A/-, 

S-     square  of  geodesic  distance  between  the  two  points  on  the  spheroid, 
r}~     pr  at  the  point  .v,  _)',, 

rnl-  -  pi'  at  the  mean  latitude  of  the  two  points  A-,  r,,  .v.,  i'.,, 
tf>  =  that  mean  latitude, 
/--•tan  6,      •//-     r'J  cos- c/>  (i  —<•-), 
A.v   .  A.vAv 


v/VA.v'r,      -//-VA.v-Ar 


OF  PROJECTIONS 

A.vAr  , 


\         ,s   a.*       A-vO'i  '^)      •'/-•/  A.v'-'  A  r 
Aa    -  o,  +  ft.,  = 


spheroidal  excess  of  space  //„///>'  7^, 


Jr  cos-  /,    -  A.v-  .,  (•?  r.-  +  6v,  Ar  +  2Ar)  •    ' 

ym  v" 

5s  si.r  /,  =  Ar  +  AA'L>,  (  ?  r,  Ar  +  A,-)  +  9lf/  ^  A-r  .. 


-  AA--  + 

A   -      A'V" 

y,,,~ 

Ay-~i( 

»  J'i~  +  3J'i  Ar  +  A>'2)  '•  -  -  » 

,                   >s       7/J/A.r'A7 

'!(^-Av 
(j'i  ^.l'J» 

Ar) 

o    r  L 

A.r  3 

i3       7/-7A.v:!Av 

which  agree  with  Helmerfs  fonnulae  remembering  tb.at  he  reckons  .v  and_>' 
as  positive  S  and  W. 

The  above  expressions,  as  well  as  those  of  Helmert,  Jordan  and 
Schreiber,  refer  to  the  gecjdesic  lines  on  the  spheroid  ami  iheir  directions. 
Before  using  them  the  appropriate  corrections  should  'he  made  to  the 
observed  directions,  or  we  may  add  these  corrections  to  the  above  expressions, 
thus  giving  the  observed  directions.  This  amounts  to  adding  ^r/-'A.vAr  /•„,- 
and  l-if  A.v  A  y'r,,,-  to  the  expressions  for  ^  and  &_,,  i.e.  instead  of  writing  Ar 
in  the  first  term  of  each  we  write  Ar  (i  -i-?/-).  As  these  corrections  are  so 
small  as  to  be  generally  neglected  we  may  well  reject  the  higher  terms  in  //-, 
and  use  only  the  terms  involving  the  second  powers  of  the  coordinates,  with 
the  proper  value  of/',,,. 

In  the  case  of  the  .Minimum  Error  (  "assim  Projection,  our  A.v  is  reduced 
by  ]  V~  '',„'  A.  v,  where  V  is  the  limiting  value  of  the  }r  coordinates.  The 
effect  of  this  on  the  bearing  of  the  direct  line  from  .v,_r,  to  ,v._,_r,,  is  to  deduct 

A.vArl'"  ,-        i 

\\e  shall  therefore  have 

A.v(2V]  4-  r.,)      A.vAr 

ft,  --  -.,   -     +  ,     .,',;  (  i'f  -  r,r-.  +  >'a  -  ]    )+••-, 

o  >'..„-  (v,,'~*~ 

A.v  (  r,  -r  jr.)      A.vA 


,     .    ,  (  )',-  -  r,  r-2+  i'-"- 
or,,,-s- 


Similarly  we  shall  have 

A 


A.v- 


\l-6r 
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If,  instead  of  making  the  average  scale  error  a  minimum,  we  make  the 
maximum  scale  error  at  the  extreme  meridian  Y  equal  and  opposite  to  the 
scale  error  along  the  principal  meridian,  as  in  the  International  Map 
on  the  scale  of  i-Million,  we  have  only  to  replace  the  Y"  in  the  above 
formulae  by  :j  Y-  to  obtain  the  corresponding  finite  errors  of  bearing 
and  distance. 

The  expressions  can  also  be  written  in  the  form 

Or  +  AJ'a  +}'•?  -  y2)  sin  /!  cos  A  ...  ], 
)  -  (yi2  +}'iy-2  +>'•/  -  Y-)  sin  A  cos  /i  .  .  .  }, 


=  s      - 


For  the  orthomorphic  projections  they  reduce  to 

Sj  =    —  (2^!  +j'a)  +  terms  in  4th  powers  of  coordinates, 


remembering  that  8]  and  &,  will  be  in  circular  measure,  and  that  if  we  want 
the  scale  error  along  the  principal  meridian  equal  and  opposite  to  that  along 
the  extreme  meridian  we  must  write  It  Y~  instead  of  Y'-. 

Piz/etti  on  p.  391  gives  formulae  for  the  error  in  distance  and  bearing 
of  a  line,  on  the  Direct  Conical  Orthomorphic  Projection,  starting  from  the 
mid-parallel  with  azimuth  «,,,  which  appear  to  be  erroneous.  For  S  -  s  he 

S"  COS"  a,,  s*  COS-  a0 

gives    -       --     which  should  read        —-r-    to  agree  with  the  above  (putting 


s'2  COS'  a,, 

v,  -  o  and  v..  =  x.>  in  direct  system  r:  s  cos  /]).      1'or  6,  he  gives  e  =  -       -— 

up,  i  /Vn 

r  sin  a,,  cos  «„ 
which  should  read  -    . 


An  error  in  Jordan's  Handbuch. 

In  checking  the  formulae  above  the  writer  had  occasion  to  use  the 
formulae  for  Reversion  of  Series  given  by  Jordan  on  p.  191,  and  in  doing 
so  found  an  error  in  the  third  line  of  terms  for  the  value  of  ,v  printed  thus: 

'•'•(f        A\f       A"' if  '    A*  if       A~'(f       A" a"'      A4 if      A:'if/ 

T.CD 

The  coefficient  of  the  fourth  term  is  in  error.      It  should  be          .,  instead 

,  4  CD 
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CHAPTER  VI 

THE  CONVERGENCY  OF   MERIDIAN'S 

The  question  of  Convergency  of  Meridians  has  always  presented  some 
difficulty  to  the  writer,  and  probably  to  others  also,  especially  when  applied 
to  the  spheroid.  Jordan  gives  two  or  three  definitions  of  it  in  different  parts 
of  Vol.  in  of  his  Handbuch,  and  yet  another  definition  at  p.  501  derived 
from  Gauss  when  dealing  with  the  Conformal  Projection.  This  latter  he- 
distinguishes  by  calling  it  the  Convergency  of  Map  Meridians,  being  the 
angle  between  the  tangent  at  a  point  of  a  projected  meridian  of  the  spheroid 
and  the  principal  meridian  of  the  map.  This  of  course  is  the  convergency 
which  a  user  of  the  map  or  projection  wants  to  know  when  converting 
a  bearing  calculated  from  the  coordinates  into  a  true  bearing.  It  will  be 
noted  however  that  the  convergency  so  derived  differs  from  the  convergency 
as  derived  from  the  Cassini  Projection  by  a  small  quantity  =  HX"?/-'  sin  4>  cos"  4>- 
An  explanation  of  this  difference  is  given  by  Jordan  on  pp.  503-4  which 
for  some  time  puzzled  the  writer  and  it  was  not  until  he  had  studied  the 
work  of  J.  I.  Craig  on  the  theory  of  the  indicatrix  that  the  explanation 
became  clear. 

It  might  appear  that  the  A-  coordinate  in  the  Gauss  Conformal  Projection 
was  the  same  as  in  the  Cassini.  It  is  so  in  the  case  of  the  sphere,  but  not 
in  the  case  of  the  spheroid.  By  expanding  the  formulae  for  the  coordinates 
on  the  Cassini  after  the  manner  of  Jordan  still  further  the  writer  finds 

.v  -  J/4-  i-rX-  sin  (f>  cos  cf>  +  .,\  i'X4  sin  </>  cos:!  <£  (5  -  t~  +  517")  +  . . . , 

y  =  i'X  cos  <b  -  }.vX':  sin'2  </>  cos  <£  —  --,  *  0r/V'  sin-  4>  cos:;  </>  (S  -  t~  +  Si/-';  —  ... . 
whereas  on  the  Gauss  Conformal  we  have  (Jordan,  p.  497) 

x  =  M+  ^  vX-  sin  </>  cos  <£  +  J,  i'X4  sin  <£  cos:;  <£  ( 5  -  t~  +  977-)  +  . .. , 

y  •=  v\  cos  <f>  +  ,V  )^X:;  cos:i  </>  ( i  —  /'-'  +  i?;) 

-  jwo-vX'"'  cos:'  (f)  (5  -  i8/-  +  /4+  1477-  —  5S-//-/-)  — 

It  will  be  observed  that  in  the  term  in  X4  for  x  we  have  +  5,;-  in  the  Cassini 
as  against  +  977-  in  the  Gauss. 

The  following  value  for  the  convergency  on  the  Cassini  is  also  found  by 
continuing  the  method  of  Jordan  : 
f       X'-cos-  <f>  . 


Now  the  work  of  J.  I.  Craig,  p.  34,  enables  us  to  find  the  true  bearing 
of  the  line  on  the  spheroid  to  which  a  line,  through  any  poini  on  the 
projection,  parallel  to  the  axis  of  .v  corresponds.  It  is 

r  cos  (/>  (/'.'  </<t> 


tan  a  -  -  -  , 

i>        ay 
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Performing  these  differentiations  we  find 
dyjd<}>  =  —  p\  sin  <£  j  i  +  -  -  (2  —  /-  +  2-rf) 


24 
whence 


.      .  f          X2COS2d>  . 

a  =  A  sin  <£  J  i  +  i  +  /-  +  r- 


X'cos4*/,,  , 

+  (2  +  At'  +  2T  +  AY)  •••  K 

15 

X4  cos4  <£  , 


agreeing  with  the  convergency. 

When  however  we  apply  the  same  method  to  the  Gauss  Conformal 
Projection  to  find  the  true  bearing  on  the  spheroid  to  which  a  line  of  zero 
bearing  on  the  projection  corresponds  we  shall  find 

-\  A      \  ^'    C°S'2    ^    /  °N  ^    C°S4    ^  r-       ^       "\  1 

Is  :s  J  ' 

differing  from  the  value  in  the  Cassini  by  the  small  quantity 
r!  rj'2  X:1  sin  <f>  cos'J  c/>  +  .... 

This  shews  that  in  the  Conformal  System  the  lines  of  zero  bearing  do  not 
correspond  with  the  lines  of  zero  bearing  in  the  Cassini  System,  or  with  the 
lines  of  same  bearing  on  the  spheroid.  The  indicatrix,  though  a  circle  in 
the  Conformal  System,  is  at  every  point  turned  through  the  above  small 
angle.  This  is  a  point  not  of  practical  importance  perhaps,  as  the  quantity 
in  question  is  so  small,  but  it  has  some  theoretical  interest  in  understanding 
the  true  nature  of  the  convergency  of  map  meridians. 

For  the  scale  value  along  a  line  parallel  to  the  principal  meridian  on 
the  projection  we  also  find  from  J.  I.  Craig's  formulae 

'dx  dy      dx  dy  \  ,//  .,  /v/v\-  fdy 


;/x    ,/A  +      p 


which  in  the  case  of  the  Cassini  leads  to 


whereas  in  the  (lauss  Conformal  we  have 

ws=.()  =  i  +  o  X-  c:os';  6  ( i  -  7/- )  +  .J1,  X;  cos'  </>  (5  -  4/-  -  i4'//-'  -  2Sv/-/-)  -f  — 
For  any  other  value  of  ft  we  have  approximately  on  the  Cassini 

111  I     ;     ,!,  X-  COS-  (J3  (  I    +  '(/-')  COS'  ft  -r   .... 
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Some  small  errors  in   Helmert's   Formulae 
for  Convergency,   etc. 

In  checking  the  formulae  for  the  convergency  on  the  spheroid,  etc.  by 
Helmert's  values  the  writer  was  led  to  discover  some  small  errors  in  the 
latter,  which  should  be  pointed  out,  though  they  will  hardly  have  any 
appreciable  effect  in  practice. 

It  will  be  found  that  the  value  for  the  convergency  can  be  put  into 
the  form 

tan  y  -  tan  A  sin  <£  j  i  +  I  if  cos-  c£A'J  +  ,'.  rj-  cos1  $A4  (  7      5/-')  +    .  .  j  , 
Transferring  this  into  Helmert's  form,  we  write 

c'-'  cos'J  Z>'., 
y  =  t,     A  ^  Z,,  ._,,      <J>  ^  />,,     -t/-  -  .,  -  , 

whence  we  have 


tan/-  tan  Z.^sinZ',  |i  +  ?.£-2Z-  cos4/*..  [i 
as  against 

tan  t     tan  Z,_  ,,  sin  /?,  \  i  +  \  f  L'-  cos4  Z',  [i  +  .V2  (5  -  3  cos  /?,) 

*•  --^z(-ff-§cos'-^)]J, 

as  given  by  Hel  inert,  p.  446,  equation  (7). 

The  error  seems  to  arise  in  the  equations  (2),  (4),  (5)  and  (6)  on  the 
previous  page  445.     In  equation  (2)  the  fourth  line  should  read 

41.rZ4  (30  —  25  shr  /'cos'-'  F—  16  cos4  F), 
instead  of  ,'.  L  '  (30  —  10  sin-  /''cos-  F—  cos4  F). 

Equations  (4)  and  (5)  are  in  error  because  the  terms  in  r  are  omitted, 
though  they  are  correctly  given  in  equation  (10),  p.  442,  for  the  relation 
between  cos  /'"and  cos/>V  These  omissions  affect  equation  (6)  and  the 
following  equation  which  should  read 

sin  /i.,  =  sin  J3.,\i  —  .1  r  cos'-  f\_i  +  -\e-  (1-3  sin-  F)  -f  L-  sin-  F 

+  -|rZ2sina/'Xi  +  sin-^)-  }.  Z4siir  F~  Z4sin'  F 
4-  ^  t'4  (i  +  2  sin-  Z''+  5  sin4  Zy")  ][•  , 

tan  /  =  -  tan  Z  sin  B  .,  \  \  +  \  f-  Z-  cos4  F"\_i  +  r  +  L-  (;",  -  |:  cos'J  /•')];  . 
Helmert's  investigation  is  extremely  difficult  to  follow  and  the  use  of  />',. 
the  actual  latitude  of  a  point  and  F  the  latitude  of  the  point  where  the 
perpendicular  coordinate  of  the  first  point  cuts  the  principal  meridian,  in 
the  same  formula,  is  very  apt  to  lead  to  confusion. 

On   the   sphere   the    latitude   of  the    foot   of    the   perpendicular    which 
determines  the  .v  coordinate  in  the  Cassini  Projection  is  given  by 

tan  F  —  tan  />'._,  sec  Z 

in  Helmert's  notation,  or  tan  <£,  -  tan  6  sec  A  in  ours.      From  this  tlu  well- 
known  expansion  for  <£,  -  <t>  is  derived  by  trigonometry. 
the  same  formula  can  be  put  in  the  form 

tan  6,  ---  tan  d>  sec  A  [  i   :  !".  //-'A'-'  (  i  -:  :';  A':  cos  26  )  .  ..  |  . 


76  EXPANSIONS 

In  the  same  way  the  formula  for  the  y  coordinate  can  be  put  in  the  form 
sin  (y/v.)  -  sin  A.  cos  </>  (i  -  T\  ?fA.4  sin2  <J>  cos2  <£). 


List  of  Expansions  used  in  these  Papers. 

sinfl^-^  +  T-^-^V^'  +  Winnr*9--", 

cos0=i-^2  +  ^-Ti^6  +  To^o08----, 
tan  0  =  6  +  -J08  +  T2T0"  +  jtfj  67  +  a-f^fl9  +  ...  , 


sec   = 


tan-1  0  =  0-  ^:f  +  i&5-  i-07  +  J  6!l-  ..., 
log,  sin  0  =  log  6  -  ^  -  Ti^  -  ^V**8  -  T3T^o  o  ^8 
log,  sec  ^  =  J-  &'  +   *s  6*  +  TV  (?  +  ^i  A  -ff  68  +  .  .  .  , 
loge  tan  ^  =  log  6  +  ^  +  ^P+^V«  +  TOT^+  -, 
\oge(i+0)  =  6-±6-  +  ±8:i-iOi  +  lF-..., 
\ogc(i-6)  =  -  (e  +  },6-  +  1&  +  \P  +  U5+  ...j  , 


tan  (6  +  //)  =  tan  0  +  h  (  \  +  /2)  +  /izf(i+  /2)  +  ^//3  (  i  +  4/2  +  3/4  )+..., 
where  /  =  tan  0. 
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